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A  one-dimensional  (1-D)  ecological  model,  HABSIM,  examined  the  initiation  and  maintenance  of  the 
2001  red  tide  on  the  West  Florida  shelf  (WFS).  Phytoplankton  competition  among  toxic  dinoflagellates 
(Karenia  brevis),  nitrogen  fixing  cyanophytes  ( Trichodesmium  erythraeum ),  large  siliceous  phytoplankton 
(diatoms),  and  small  non-siliceous  phytoplankton  (microflagellates)  explored  the  sequence  of  events 
required  to  support  the  observed  red  tide  from  August  to  December  2001.  The  ecological  model 
contained  24  state  variables  within  five  submodels:  circulation,  atmospheric  (iron  deposition),  bio¬ 
optics,  pelagic  (phytoplankton,  nutrients,  bacteria,  zooplankton,  and  fish),  and  benthic  (nutrient 
regeneration).  The  2001  model  results  reaffirmed  that  diazotrophs  are  the  basis  for  initiation  of  red 
tides  of  K.  brevis  on  the  WFS.  A  combination  of  selective  grazing  pressure,  iron  fertilization,  low  molar 
nitrogen  to  phosphorus  ratios,  and  eventual  silica  limitation  of  fast-growing  diatoms  set  the  stage  for 
dominance  of  nitrogen  fixers.  “New”  nitrogen  was  made  available  for  subsequent  blooms  of  K.  brevis 
through  the  release  of  ammonium  and  urea  during  nitrogen  fixation,  as  well  as  during  cell  lysis,  by  the 
Trichodesmium  population.  Once  K.  brevis  biomass  reached  ichthyotoxic  levels,  rapid  decay  of 
subsequent  fish  kills  supplied  additional  organic  nutrients  for  utilization  by  these  opportunistic  toxic 
algae.  Both  nutrient  vectors  represented  organic  non-siliceous  sources  of  nitrogen  and  phosphorus, 
further  exacerbating  silica  limitation  of  the  diatom  population.  The  model  reproduced  this  spring 
transition  from  a  simple  estuarine-driven,  diatom-based  food  chain  to  a  complex  summer-fall  system 
of  Trichodesmium  and  toxic  dinoflagellates.  While  the  model  was  able  to  replicate  the  initiation  and 
maintenance  of  the  2001  red  tide,  bloom  termination  was  not  captured  by  this  1-D  form  on  the  WFS. 
Here,  horizontal  advection  and  perhaps  cell  lysis  loss  terms  might  play  a  significant  role,  to  be 
addressed  in  future  three-dimensional  simulations. 

©  2012  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Harmful  algal  blooms  (HAB)s  of  the  toxic  dinoflagellate, 
Karenia  brevis,  occur  annually  along  the  West  Florida  coast 
(Steidinger  et  al„  1998)  with  far-held  implications  for  both  the 
South  Atlantic  Bight  and  northern  Gulf  of  Mexico  (Walsh  et  al., 
2006).  Origins  of  these  red  tide  blooms  have  been  traced  to 
bottom  initiation  along  the  20-30  m  isobaths  on  the  West  Florida 
shelf  (WFS)  between  Cedar  Key  and  Naples  (Steidinger  et  al., 
1998;  Walsh  et  al.,  2009;  Weisberg  et  al.,  2009)  in  the  eastern  Gulf 
of  Mexico  (COM).  Subsequent  physical  transport  within  the 
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bottom  Ekman  layer  (Weisberg  et  al.,  2009)  plays  a  critical  role 
in  the  timing  and  location  of  the  blooms  relative  to  both  nutrient 
sources  and  grazing  stresses  (Walsh  et  al.,  2009,  2011a).  Hence, 
the  economic  and  health  impact  of  local  red  tides  to  adjacent 
human  population  centers  depends  greatly  upon  both  the  avail¬ 
ability  of  nutrients  and  herbivores  during  physical  transport  of 
the  bloom. 

Based  upon  recent  data  from  175  cruises  on  the  WFS  during 
March  1998-December  2001  (Walsh  and  Steidinger,  2004),  we 
arrived  at  a  consensus  on  a  repetitive  sequence  of  required  events 
for  annual  red  tide  formation  in  the  GOM  (Walsh  et  al.,  2006).  Weak 
upwelling  by  the  Loop  Current  at  the  shelf-break  (Weisberg  and  He, 
2003)  provides  Caribbean  seed  populations  of  the  N2  fixing  cyano- 
phyte,  Trichodesmium  erythraeum,  within  a  nutrient  milieu  of  low 
nitrogen/phosphorus  (N/P)  ratios  on  the  WFS.  Background  K.  brevis 
populations  of  ~0.01  pg  chi  1_1  occur  throughout  the  GOM  (Geesey 
and  Tester,  1993).  The  low  N/P  ratios  of  <  6  on  the  central  WFS  are 
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the  consequence  of  (1)  estuarine  supplies  of  leached  Miocene 
phosphorite  deposits,  and  (2)  different  rates  of  element  recycling 
of  plankton  debris  (Walsh  et  al.,  2006). 

Within  the  WFS,  these  phosphorus-rich  nutrient  supplies 
modulate  seasonal  phytoplankton  succession  after  the  initial 
diatom  spring  bloom,  once  summer  deposition  events  of  iron¬ 
laden  Saharan  dust  allow  Trichodesmium  blooms  to  utilize  ubi¬ 
quitous  dissolved  dinitrogen  gas  (N2)  of  otherwise  nitrogen-poor 
seawater  near  the  shelf-break  (Lenes  et  al.,  2001;  Walsh  and 
Steidinger,  2001).  The  ambient  low  stocks  of  nitrate  (Masserini 
and  Fanning,  2000)  and  selective  grazing  stresses  exerted  by 
copepods  (Kleppel  et  al.,  1996)  place  the  diatoms  at  a  competitive 
disadvantage  here,  compared  to  the  slower  growing  diazotrophs 
and  K.  brevis  (Walsh  et  al.,  2001a,  2003,  2006).  These  co-occurring 
non-siliceous  phytoplankton  populations  initially  reside  in  the 
bottom  Ekman  layer  due  to  similar  diel  vertical  migration 
patterns,  with  the  sun-adapted  diazotrophs  seeking  higher  phos¬ 
phate  and  the  shade-adapted  dinoflagellates  seeking  lower  light 
(Walsh  et  al.,  2006).  Fall  upwelling  events  both  resurface  these 
phytoplankton  residents  of  the  bottom  Ekman  layer  higher  into 
the  euphotic  zone  and  transport  them  shoreward  (Walsh  et  al„ 
2003;  Weisberg  and  He,  2003;  Weisberg  et  al.,  2009),  such  that 
shelf-shading  components  of  these  toxic  dinoflagellate  popula¬ 
tions  can  be  finally  detected  by  satellites. 

In  contrast,  the  faster-growing  diatoms  otherwise  win  the 
spring  phytoplankton  competition  within  COM  coastal  waters 
(Gilbes  et  al.,  1996;  McPherson  et  al„  1996;  Phlips  and  Badylak, 
1996;  Dortch  et  al.,  1997;  Bledsoe  and  Phlips,  2000;  Gilbes  et  al„ 
2002;  Qian  et  al.,  2003;  Jurado  et  al.,  2007;  Heil  et  al.,  2007), 
where  both  dissolved  nitrate,  silicate,  and  iron  stocks,  as  well  as 
herbivores,  are  sometimes  abundant  in  the  estuaries,  until  sea¬ 
sonal  silica-limitation  ensues  here  and  on  most  shelves  (Walsh 
et  al.,  2011b).  Farther  offshore,  the  usual  nitrate/silicate  ratios  do 
not  favor  diatom  competitors  on  the  WFS  (Walsh  et  al.,  2003). 
Moreover,  both  initial  excretion  and  then  lysis  of  the  offshore 
Trichodesmium  populations  during  their  bloom  decay  can  provide 
a  unique  nutrient  source  for  red  tides  on  the  outer  shelf  (Walsh 
and  Steidinger,  2001;  Lenes  et  al.,  2001,  2008;  Mulholland  et  al„ 
2004,  2006;  Walsh  et  al.,  2006;  Lenes  and  Heil,  2010),  which  the 
non-migratory  diatoms  evidently  do  not  access. 

These  recycled  nutrients  in  the  form  of  dissolved  organic  nitro¬ 
gen  (DON),  i.e.,  urea  favored  by  K.  brevis  (Steidinger  et  al.,  1998),  and 
dissolved  organic  phosphorus  (DOP),  can  be  utilized  by  these  toxic 
dinoflagellates  (Vargo  and  Howard-Shamblott,  1990),  giving  them  a 
competitive  advantage  over  faster  growing  phytoplankton  that 
prefer  inorganic  nutrient  sources.  Once  seasonal  K.  brevis  biomass 
has  exceeded  ~1  x  105  cells  1  i.e.,  concentrations  great  enough  to 

cause  fish  kills  (Landsberg,  2002),  high  fish  decomposition  rates  of 
up  to  50%  per  day  prevail  within  warm  summer  Florida  waters 
(Stevenson  and  Childers,  2004),  quickly  converting  particulate 
nutrients  of  fish  flesh  to  the  dissolved  pools. 

Walsh  et  al.  (2006)  estimated  that  there  was  ~  1  individual  fish 
per  m2  as  photographed  within  the  August  1947  drift  lines  of 
coincident  dead  fish,  diazotrophs,  and  toxic  dinoflagellates  (Gunter 
et  al.,  1948).  These  fish  populations  were  dominated  by  white 
grunts  ( Haemulon  plumieri),  pinfish  (Lagodon  rhomboides ),  and 
thread  herring  ( Opisthonema  oglinum).  Previous  estimates  of  mean 
fish  wet  weight  (ww)  of  90  g,  dry  weight  (dw)  of  25%  ww, 
N  content  of  ~10%  dw,  and  P  content  of  0.7%  dw  suggested  a 
potential  nutrient  supply  for  the  maintenance  of  a  red  tide  of 
~16  pmol  PN  kg-1  and  ~0.5  pmol  PP  kg-1  within  a  10-m  eupho¬ 
tic  zone  (Walsh  et  al.,  2006). 

In  addition  to  decomposing  cyanophytes  and  dead  fishes, 
several  other  nutrient  sources  such  as  benthic  flux  (Darrow 
et  al.,  2003),  ground  water  loading  (Hu  et  al.,  2006),  riverine 
discharge  (Brand  and  Compton,  2007),  and  zooplankton  excretion 


(Lester,  2005;  Lester  et  al.,  2008)  may  all  contribute  to  the 
development  and  maintenance  of  K.  brevis  blooms.  While  the 
relative  role  of  each  nutrient  source  is  still  in  question,  it  is  clear 
that  K.  brevis  is  an  opportunistic  organism  capable  of  utilizing  all 
available  resources  (Vargo  et  al.,  2008). 

Therefore,  continued  development  of  validated  simulation 
models  of  those  environmental  factors  which  lead  to  interannual 
successes  and  failures  of  “red  tides”  of  K.  brevis  becomes  increas¬ 
ingly  important  in  eventual  operational  forecasts  of  their  spatio- 
temporal  bloom  dynamics  (Walsh  et  al.,  2011c),  by  teasing  out 
the  relative  role  of  each  nutrient  source  in  a  heuristic  mode.  Of 
greater  importance,  however,  will  be  eventual  prudent  manage¬ 
ment  strategies  for  this  and  other  pending  future  HABs  of  Florida 
waters  (Walsh  et  al.,  2011a). 

Over  the  last  50  years,  at  least  15  simple  (Slobodkin,  1953)  and 
increasingly  more  complex  (Walsh  et  al.,  2001a,  2002,  2003;  Jolliff, 
2004;  Milroy  et  al.,  2008;  Lenes  et  al.,  2008)  models  of  red  tides  in 
the  eastern  COM  have  explored  these  earlier  hypotheses,  with  some 
numerical  formulations  relying  on  synoptic  estimates  of  satellite 
color  signals  for  validation  data  (Jolliff  et  al.,  2003;  Walsh  et  al., 
2006,  2009).  We  also  used  this  remotely  sensed  surface  information 
for  initial  and  boundary  conditions,  providing  constraints  on  which 
model  parameters  can  reproduce  both  the  seasonal  succession  of 
biochemically  mediated  dominant  groups  of  phytoplankton  and  the 
complex  ecological  trophic  cascade  for  initiation  of  more  recent 
larger  red  tides  of  (1)  K.  brevis  and  (2)  a  competing,  previously 
cryptic  saxitoxic  dinoflagellate  Pyrodinium  bahamense  in  the  GOM 
(Walsh  et  al.,  2006,  2011  a). 

People-killing  HABs  of  the  latter  P.  bahamense  were  involved  in 
more  than  100  deaths  of  humans,  upon  consumption  of  saxitoxin- 
adulturated  sea  food  from  the  South  China  and  eastern  Caribbean 
Seas  over  the  last  few  decades  (Walsh  et  al.,  2011a).  Future 
simulation  studies  must  thus  address  what  are  the  “rules  of 
engagement”  for  relatively  benign  ichthyotoxic  K.  brevis  on 
the  WFS  to  be  replaced  there  and  downstream  by  saxitoxic 
P.  bahamense,  whose  poisons  can  be  much  more  deadly  to  humans 
than  sarin,  curare,  and  brevetoxins. 

Yet,  the  simpler  objective  of  the  present  study  was  to  evaluate 
the  fidelity  of  an  already  complex,  multi-species  biochemical 
model  for  simulation  of  those  ecological  processes,  which  were 
required  for  the  initiation  and  maintenance  of  K.  brevis  HABs  on 
the  WFS  (Fig.  1).  The  model  was  thus  quasi-independent  of  the 
three-dimensional  (3-D)  circulation  and  ignored  any  dinoflagel¬ 
late  competition  between  brevetoxic  and  saxitoxic  poison 
sources.  Therefore,  we  present  the  hindcasts  for  the  2001  red 
tide  of  K.  brevis,  utilizing  the  one-dimensional  (1-D)  form  of  the 
Harmful  Algal  Bloom  SIMulations  (HABSIM)  model  in  the  vertical 
dimension  to  test  its  fidelity  with  four  different  cases,  calibrated 
against  prior  in  situ  cell  counts  and  subject  to  a  sensitivity 
analysis  of  the  model’s  parameters.  We  then  discuss  needed 
additions  to  bring  to  fruition  operational  forecasts  of  WFS  HABs. 

Accordingly,  even  more  complex  future  3-D  HAB  forecasts 
would  then  involve  revised  ecological  formulations  of  time- 
dependent  selective  nutrient  regeneration,  based  on  the  present 
simulation  analysis.  But,  to  predict  the  outcome  of  competition  of 
the  relatively  benign  fish-killer  K.  brevis  against  the  very  poiso¬ 
nous  P.  bahamense  to  humans  under  nutrient-depleted  conditions 
(Walsh  et  al.,  2011a),  spatial  teleconnections  of  these  two  dino¬ 
flagellate  HABs  must  also  be  explicitly  considered  within  the  shelf 
waters  and  estuaries  of  the  GOM  and  the  downstream  South 
Atlantic  Bight. 

Future  answers  to  complex  questions  about  the  possible 
impacts  of  legislated  changes  of  fertilizer  applications  to  lawns 
and  agricultural  crops  unfortunately  require  complex,  coupled 
biophysical  models.  Here,  we  evaluate  the  performance  of  HABSIM 
with  95  parameters,  providing  insight  to  the  requirements  of  the 
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Fig.  1.  Location  of  the  one-dimensional  model  (X)  west  of  Sarasota  above  the  20-m  isobath  in  relation  to  Tampa  Bay  and  Charlotte  Harbor.  The  bounded  area  represents 
the  region  over  which  the  maximum  weekly  Karenia  brevis  concentrations  were  obtained  during  2001 ,  for  validation  of  the  present  model. 


final  generation  of  operational  HAB  models  of  the  southeastern 
United  States,  within  a  forthcoming  era  of  ecosystem-based  man¬ 
agement  of  marine  biotic  resources,  now  perturbed  by  overfishing, 
rising  water  temperatures,  and  oil  spills. 

2.  Methods 

The  present  HABS1M  version  consisted  (Fig.  2)  of  24  state 
variables:  temperature,  salinity,  kz  for  vertical  mixing,  spectral 
light,  colored  dissolved  organic  matter  (CDOM),  dissolved  organic 
carbon  (DOC),  dissolved  organic  nitrogen  (DON),  dissolved 
organic  phosphorus  (DOP),  dissolved  inorganic  carbon  (DIC), 
nitrate + nitrite  (N03),  ammonium  (NH4),  phosphate  (P04),  silicate 
(Si04),  iron  (Fe),  four  functional  groups  of  phytoplankton  (dia¬ 
toms,  microflagellates,  T.  erythraeum,  and  K.  brevis ),  fecal  pellets, 
ammonifying  and  nitrifying  heterotrophic  bacteria,  detritus, 
larval  fish,  adult  fish,  and  sediment  microbiota. 

The  model  had  five  submodels  of  (1)  circulation,  (2)  dust 
deposition,  (3)  bio-optics,  (4)  pelagic  ecological  processes,  and  (5) 
benthic  ecological  processes,  of  which  the  last  four  compose 
HABS1M  (Fig.  2).  The  last  two  were  forced  by  outputs  of  the  first 
three  submodels.  HABS1M  was  built  upon  the  structure  of  pre¬ 
vious  ecological  models  (Walsh  et  al.,  2001  a, b,  2002,  2003,  2006; 
Jolliff  et  al.,  2003;  Jolliff,  2004;  Lenes  et  al.,  2008;  Darrow,  2008; 
Milroy  et  al.,  2008).  Therefore,  we  limit  the  methods  section  to 
description  of  the  major  modifications  and  additions.  A  complete 
set  of  equations  is  presented  in  Section  6,  Math  formulae,  while 
the  model  parameters  are  in  Table  1 . 

2.2.  Physical  processes 

The  model  had  a  1-m  vertical  resolution  over  20  depth  levels  of 
the  water  column,  with  one  sediment  layer  of  1  cm  thickness.  A  30  s 
time  step  was  used,  beginning  on  1  January  2001.  The  model  was 


spun-up  for  3  months  with  simulation  results  shown  from  April- 
December  2001.  It  was  positioned  at  the  20-m  isobath  (ECOHAB 
station  29  at  ~  27.1 7  N,  82.9CTW)  southwest  of  Sarasota,  FL  (Fig.  1). 
Physical  mixing  of  the  model  in  the  vertical  direction,  kz,  was 
calculated  with  a  Mellor-Yamada  Level  2.5  turbulence  closure 
scheme  (Mellor  and  Yamada,  1982),  in  which  daily  average  wind 
speeds  (Fig.  3(a))  were  obtained  from  the  National  Data  Buoy  Center 
(NDBC)  observations  ~195  km  WNW  of  Tampa,  FL. 

Daily  temperature  and  salinity  profiles  (Fig.  3(b)  and  (c))  were 
calculated  for  the  study  site  (Fig.  1 )  from  the  WFS  circulation  model 
(Barth  et  al.,  2008;  Weisberg  et  al.,  2009),  which  is  a  regional 
application  of  ROMS  (Regional  Oceanic  Modeling  System,  e.g., 
Shchepetkin  and  McWilliams,  2005)  nested  in  HYCOM  (HYbrid 
Coordinate  Ocean  Model,  e.g.,  Chassignet  et  al.,  2003).  As  such,  it  is 
fully  3-D  and  baroclinic,  driven  by  complete  physical  forcing  fields 
(deep-ocean  fluxes  across  the  open  boundary  from  HYCOM  and 
local  momentum  and  buoyancy  fluxes  over  the  eastern  COM  and 
WFS)  with  approximate  2-5  km  horizontal  resolution. 

The  internal  light  fields,  as  well  as  the  water-leaving  radiances, 
were  computed  from  the  bio-optical  submodel’s  inherent  optical 
properties.  They  were  a  function  of  some  of  the  absorption  and 
backscatter  properties  of  seawater  and  the  state  variables  of  the 
ecological  model,  i.e.,  total  phytoplankton  biomass,  total  hetero¬ 
trophic  bacteria,  CDOM,  particulate  detritus,  and  inorganic  suspended 
sediments  (Jolliff,  2004).  Light  at  any  given  depth  was  determined 
using  the  Beer-Lambert  law  of  exponential  decay,  with  appropriate 
attenuation  coefficients  for  blue  and  red  light  determined  as  a 
function  of  the  concentrations/stocks  of  CDOM  and  equivalent 
chlorophyll  biomass  of  the  phytoplankton  within  the  overlying  water. 

2.2.  Biological  processes 

2.2.1.  Phytoplankton  competition 

Debate  over  the  feasibility  of  multi-plankton  functional  types 
(PFTs)  and  enhanced  complexity,  versus  simpler  model  formulations, 
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Fig.  2.  Schematic  of  the  HABSIM  model  with  a  color  coded  breakdown  of  the  five  components:  green  — pelagic  phytoplankton;  orange— physical  model;  yellow=optical 
model;  purple = nutrient  pathways;  and  red  — loss  to  higher  trophic  levels.  The  blue  area  represents  the  water  column  and  the  gray  are  represents  the  benthos.  (For 
interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


continues  to  rage  within  the  oceanographic  community  (Anderson, 
2005;  Flynn,  2005;  Hood  et  al„  2006;  Friedrichs  et  al.,  2007,  2009). 
While  all  choices  have  trade-offs,  it  is  critical  that  model  variables  and 
parameterizations  are  built  upon  the  insights  of  simpler  models 
(Thingstad  et  al„  2010).  Given  the  inadequate  results  of  a  decade  of 
such  WFS  models,  we  added  a  higher  level  of  complexity  in  conjunc¬ 
tion  with  the  four  PFTs  to  capture  the  annual  planktonic  succession 
and  multiple  trophic  pathways  that  lead  to  the  dominant  HABs  of 
toxic  dinoflagellates  on  the  WFS  (Walsh  et  al.,  2011a). 

On  most  continental  shelves,  microflagellates  follow  the  spring 
diatom  bloom  within  a  typical  seasonal  succession  of  phytoplank¬ 
ton,  as  functions  of  light  availability,  water  column  stability, 
nutrient  resources,  gravitational  settling,  and  herbivore  grazing 
demands.  Furthermore,  on  the  WFS,  when  episodic  strong  upwel- 
ling  at  a  N03/P04  ratio  of  ~16  from  the  deep  GOM  prevails, 
nitrate-assimilating  diatoms  again  dominate  in  both  the  “real 
world”  and  in  models  (Walsh  et  al.,  2001a,  2001b,  2003). 

Similarly,  the  usual  spring  succession  of  first  diatoms  can  be 
fueled  by  nutrient  supplies  from  northern  estuaries  in  a  N/P  ratio 
>16  (Walsh  et  al.,  2006).  For  example,  these  siliceous  phyto¬ 
plankton  win  within  March  plumes  of  the  northern  Mississippi 
(Dortch  and  Whitledge,  1992;  Fahnenstiel  et  al.,  1995),  Apalachi¬ 
cola  (Gilbes  et  al.,  1996,  2002),  Suwannee  (Bledsoe  and  Phlips, 
2000)  Rivers,  and  to  less  extent  of  the  southern  Peace  (McPherson 
et  al.,  1996),  and  Shark  (Jurado  et  al.,  2007)  Rivers  until  the 
summer  rains  begin,  over  near  shore  regions  of  the  eastern  GOM. 

Zooplankton  stocks  respond  on  the  WFS  by  maximal  accumu¬ 
lation  of  biomass  in  May,  following  a  spring  peak  of  secondary 
production  of  WFS  herbivores  associated  with  low  surface  sali¬ 
nities,  i.e„  <  36.0  found  above  the  30-m  isobath  of  the  Florida 


Middle  Ground  (FMG)  in  1969-1970  at  ~28.50GN,  84.25  W 
(Austin  and  Jones,  1974).  Another  shelf-wide  time  series  of 
zooplankton  surveys  during  1972-1974  subsequently  found  the 
same  amount  of  spring  maximal  zooplankton  biomass  of 
~250ml  displacement  volume  1000  m-3  throughout  the  water 
column  on  the  FMG,  i.e„  to  the  northwest  of  the  study  site  (Fig.  1 ), 
during  May  1970  and  May  1973,  1974  (Austin  and  Jones,  1974; 
Houde  and  Chitty,  1976;  as  redrawn  in  Walsh  et  al.,  2011c). 

This  FMG  location  is  within  the  recurrent  seasonal  chlorophyll 
plumes  of  spring  diatoms  on  the  northern  WFS,  derived  from  local 
nutrient  supplies  of  the  Suwannee  and  Apalachicola  Rivers,  as 
observed  during  March  1992  (Gilbes  et  al„  1996)  and  on  3-27 
March  2001.  A  few  months  later,  however,  an  inverse  north- 
south  gradient  is  observed  for  the  spring  zooplankton  standing 
stocks,  with  50%  more  of  the  metazoan  biomass  found  south,  than 
north,  of  Tampa  Bay  during  May  1973  and  1974  (Houde  and 
Chitty,  1976).  Such  spatial  spring  differences  in  accumulation  of 
planktonic  pi'edators  and  algal  prey  on  the  WFS  reflect  net 
population  growth  increments  of  these  herbivores,  as  well  as 
both  settling  and  grazing  losses  of  their  upstream  phytoplankton 
prey,  over  upstream  time  and  space. 

A  north-south  gradient  of  organic  matter  of  plankton  debris 
accumulated  each  year  in  the  WFS  surficial  sediments,  with  as 
much  as  1-2%  dry  weight  organic  carbon  found  near  the  FMG 
(Walsh  and  Dieterle,  1994).  But,  little  organic  matter  of  <  0.5  dry 
weight  organic  carbon  was  instead  deposited  in  sediments  of  the 
central  WFS,  since  diatom  fallout  used  to  be  only  a  major 
component  on  the  northern  shelf  (Darrow  et  al„  2003),  as  a 
consequence  of  local  freshwater  nutrient  loadings  (Gilbes  et  al., 
1996,  2002),  prevailing  circulation  patterns  (Weisberg  et  al„  1984, 
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Table  1 

The  model  variables  and  parameters. 


Description 

Value 

Units 

Pi. 2.3.4 

Diatoms,  Microflagellates,  Trichodesmium,  I<.  brevis 

mmol  C  m-3 

Ba.„ 

Ammonifying  bacteria,  Nitrifying  bacteria 

mmol  C  m-3 

Fed 

Dissolved  Iron 

mmol  Fem-3 

Fec 

Colloidal  Iron 

mmol  Fem-3 

DIC 

Dissolved  Inorganic  Carbon 

mmol  C  m-3 

DOC 

Dissolved  Organic  Carbon 

mmol  C  m-3 

DON 

Dissolved  Organic  Nitrogen 

mmol  N  m-3 

DOP 

Dissolved  Organic  Phosphorus 

mmol  P  m-3 

nh4 

Ammonium 

mmol  N  m-3 

no3 

Nitrate 

mmol  Nm'3 

P04 

Phosphate 

mmol  P  m-3 

Si04 

Dissolved  Silica 

mmol  Si  m-3 

D 

Detritus 

mmol  C  m-3 

Z 

Fecal  Pellets 

mmol  C  m-3 

F 

Fish 

mmol  C  m-3 

(C:N)i  ,2,3,4.a,n 

Nitrogen  to  Carbon  ratio  of  phytoplankton  and  bacteria 

6.63,6.63,6,6.63,5,5 

(C:P)i,2,3,4,a,n 

Phosphorus  to  Carbon  ratio  of  phytoplankton  and  bacteria 

106,106,122,106,53,53 

(C:Fe)1>2.3,4,a.n 

Iron  to  Carbon  ratio  of  phytoplankton  and  bacteria  (xlO3) 

15,25,2,20,15,15 

(C:Chl)1234 

Carbon  to  chlorophyll  ratio  of  phytoplankton 

50,100,200,30 

T 

Temperature 

°C 

Tr 

Transport 

m2  s-1 

a 

Vertical  coordinate 

Wi,3,4 

Vertical  migration/sinking  rate  for  phytoplankton 

0.5, 0.0,*,* 

m  day-1 

^max3,4 

Maximum  vertical  migration  rate  for  P3,  P4 

77,23 

m  day-1 

Sl,2,3,4,a,n 

Growth  rate  of  phytoplankton  and  bacteria 

day'1 

71.2.3.4 

Grazing  rate  on  phytoplankton 

0.02,0.06,0.02,0.01 

day'1 

£1,2,3,4.11,3^ 

Fraction  of  growth  respired 

0.1,0.15,0.5,0.25,0.7,0.7,0.3 

'I'l 

Fraction  of  diatom  growth  excreted 

0.04 

ll.  2.3,4 

Lysis  rate  of  phytoplankton 

0.008, 0.2, *,0.01 

day'1 

ma.n 

Mortality  rate  of  ammonifying  bacteria,  nitrifying  bacteria 

0.086,0.086 

day'1 

PFe.DON 

Photolysis  rate  of  colloidal  iron,  DON 

2.01,0.00017 

mmol  m-3  day'1 

a 

Scavenging  rate 

mmol  m-3  day'1 

Kfe 

Dissolved  iron  from  atmospheric  dust 

Ql,  2,3,4 

Fraction  of  lysed  phytoplankton  to  detritus  pool 

0.1, 0.1, 0.1, 0.1 

I<d 

Decay  rate  of  detritus 

0.13 

day"1 

Wd.Z 

Sinking  rate  of  detritus,  fecal  pellets 

0.1,100 

m  day-1 

Al, 2,3,4 

Fraction  of  grazed  phytoplankton  to  fecal  pellets 

0.5,0.15,0,0.5 

bl,  2,3,4 

Fraction  of  grazed  phytoplankton  to  fish 

0.075,0.045,0,0.075 

kz 

Fecal  pellet  decay  rate 

0.13 

day-1 

f 

Fraction  of  grazed  phytoplankton  excreted  by  fish  as  fecal  pellets 

0.6 

kdiss 

Dissolution  rate  of  siliceous  detritus/fecal  pellets 

0.6 

day-1 

°1. 2.3,4 

Fraction  of  lysed  material  to  DOM 

0.9, 0.9, 0.9, 0.9 

jl, 2,3,4 

Fraction  of  grazed  material  to  DOM 

0.2,0.12,0.55,0.2 

kfd 

Decay  rate  of  fish 

day-1 

kh 

coefficient  of  eddy  diffusivity 

m2s-’ 

Iz 

PAR  at  depth 

Isatl,  2,3,4 

Light  saturation  intensity  for  each  phytoplankton  group 

190,275,400,65 

pH  m-2  s_1 

fll.2,3,4,a,n 

Temperature  specific  maximum  growth  rate  for  phytoplankton 

2.2,1.8,0.7,0.8,13,1.3 

day-1 

knh4  -  l,2,4,n 

half  saturation  constant  for  NH4 

1.5,0.3,0.5,0.25 

mmol  m-3 

kno3  -  1,2,4 

half  saturation  constant  for  N03 

1.0,03,3.0 

mmol  m-3 

kpo4  -  l,2,3,4,n 

half  saturation  constant  for  P04 

0.2, 0.2, 0.1, 0.2, 0.1 

mmol  m-3 

ksio4  -  1 

half  saturation  constant  for  Si04 

0.5 

mmol  m-3 

kDON  -  4, a 

half  saturation  constant  for  DON 

1.5, 1.5 

mmol  m-3 

kDOP  -  3,4,a 

half  saturation  constant  for  DOP 

0.1.0.4.0.5 

mmol  m-3 

kFed  -  1, 2,3,4 

half  saturation  constant  for  Dissolved  Fe 

5,5,500,5 

10"6  mmol  m~3 

Omax 

maximum  iron  scavenging  rate 

0.62 

day-1 

K 

iron  scavenging  coefficient 

0.0015 

l,max 

maximum  Trichodesmium  lysis  rate 

0.3 

day-1 

K 

fraction  of  maximum  Trichodesmium  growth  due  to  nutrient  limitation 

2009),  and  prior  herbivore  abundances  (Houde  and  Chitty,  1976), 
even  before  long-term  reductions  of  grazing  stresses  as  part  of  a 
recent  trophic  cascade  (Walsh  et  al.,  2011a). 

After  sinking  of  the  ungrazed  parts  of  the  estuarine-driven 
“regular”  spring  diatom  bloom  here  (Hitchcock  et  al.,  2000), 
sequestration  of  remineralized  ammonium  by  benthic  microflora 
instead  triggered  6-fold  seasonal  increments  of  sediment  chlor¬ 
ophyll  (Darrow  et  al.,  2003).  Increases  of  the  microflagellate 
stocks  also  occurred  (Walsh  et  al.,  2003)  in  response  to  the 
recycled  nutrients,  not  captured  by  the  phytobenthos.  Since 
protozoans  heavily  grazed  the  microflagellates  on  the  WFS 


(Kleppel  et  al„  1996),  the  recycled  nutrients  were  quickly  trans¬ 
ferred  up  the  trophic  ladder. 

This  removal  process  cleared  the  phytoplankton  stage 
(Hutchinson,  1965)  for  summer  fertilization  of  poorly  grazed, 
slow  growing  diazotrophs  and  fall  red  tides  in  the  low  N/P  waters 
of  the  WFS  (Walsh  et  al.,  2003,  2009).  To  simulate  the  conse¬ 
quences  of  the  spring  nutrient  supply  from  the  northern  COM 
rivers,  a  nutrient  pulse  of  1.0  pmol  N03  kg-1,  0.1  pmol  P04  kg-1, 
and  2.0  pmol  Si04  kg-1  was  added  to  the  upper  10-m  of  the 
model’s  20-m  water  column  on  16  March  2001  of  the  baseline 
case  of  the  model. 
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Fig.  3.  The  2001  mean  daily  (A)  wind  speed  (ms  ’)  at  the  National  Data  Buoy  Center  (NDBC)  at  station  42036, 195  km  WNW  of  Tampa,  FL,  as  well  as  the  (B)  temperature 
(°C)  and  (C)  salinity  fields  calculated  from  the  WFS  ROMS  circulation  model  versus  depth  (Barth  et  al.,  2008). 


Fig.  4.  The  2001  daily  influx  of  total  dissolved  iron  (pmolm  2  day  ')  to  the  West  Florida  shelf  surface  waters  as  a  perturbation  condition.  Atmospheric  dust 
concentrations  at  Miami  and  daily  rainfall  at  Tampa  International  Airport  were  used  to  calculate  both  wet  and  dry  deposition  (after  Lenes  et  al„  2008). 


2.2.2.  Element  stoichiometry 

To  initiate  nitrogen-fixation,  iron  of  Saharan  dust  origin  was 
injected  at  the  surface  of  the  model  in  this  first  case  during  mainly 
June-July  2001  (Fig.  4).  In  each  case  of  the  model  (with  and 
without  estuarine,  diazotroph,  and  dead  fish  sources  of  nutrients), 
phytoplankton  competition  among  the  PFTs  was  effected  by 
variations  in  rate  parameters,  such  as:  the  temperature-depen¬ 
dent  maximal  growth;  photosynthesis  at  saturation  light  inten¬ 
sity;  half-saturation  constants  for  uptake  of  nutrients; 
respiration;  excretion;  grazing;  settling;  and  migration  for  each 
group  of  phytoplankton  and  benthic  microalgae  (Table  1). 

In  this  complex  ecological  model,  the  carbon-based  growth 
rates  of  pelagic  and  benthic  diatoms  required  nitrate  and  silicate 
in  a  constant  molar  ratio  of  100  C/16  N/16  Si/1  P  (Table  1)  for 
balanced  growth.  The  other  PFTs  had  the  same  N/P  Redfield  molar 
body  ratios  of  1 6/1  (Table  1 ),  such  that  the  model’s  phytoplankton 
decomposition  and  recycling  through  the  bacterial  pool  occurred 
at  this  ratio  of  10/1,  not  that  of  selective  remineralization 
of  phosphorus  (Grill  and  Richards,  1964;  Kamatani,  1969).  These 


bacterial  recycling  processes  instead  yielded  an  initial  dissolved 
N/P  ratio  of  ~4,  observed  in  benthic  depocenters  of  shallow 
waters,  downstream  of  river-grown  diatoms  on  the  northern 
GOM  (Walsh  et  al.,  2006). 

Neglect  of  this  separate  time-dependence  of  element  recycling 
of  phytoplankton  debris  in  the  model  led  eventually  to  phos¬ 
phorus-limitation  and  thus  unrealistic  time  lags  of  the  simulated 
Trichodesmium  and  subsequent  Karenia  blooms.  While  presently 
of  heuristic  interest,  such  a  failure  in  an  operational  model  would 
falsely  predict  both  late  arrivals  of  K.  brevis  FIABs  on  the  beaches 
and  inadequate  warnings  of  associated  fish  kills  and  asthma 
attacks. 

The  assumed  C/chl  weight  ratios  of  30-50  for  diatoms,  instead 
of  ~200  found  for  all  particulate  matter  during  the  NEGOM  and 
ECOHAB  cruises  (Jochens  and  Nowlin,  1999;  Heil  et  al„  2004), 
yielded  a  particulate  PN/chl  ratio  (pmol/pg)  of  ~0.6  for  siliceous 
phytoplankton,  as  measured  during  other  diatom  blooms  (Walsh 
et  al.,  1978).  Finally,  a  composite  chlorophyll  content  of 
~0.5  x  10-6  pgcell-1  implied  a  mixed  diatom  population  of  90% 
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small  (8  pm)  and  10%  large  (50  pm)  forms  of  these  netplankton 
(Walsh  et  al.,  2001a),  with  a  maximal  growth  rate  of  ~2.2  day-1  at 
30  DC  (Table  1 ).  They  had  an  intermediate  light  saturation  intensity 
of  190  pmol  photons  m~2  s_1. 

The  second  functional  group  of  the  model  was  small  microfla- 
gellates  of  6-8  pm  diameters,  with  a  C/chl  ratio  of  ~100,  and  a 
chlorophyll  content  of  ~0.5  x  10  (>  pg  cell  1  (Table  1).  They  had  a 
similar  maximal  growth  rate  of  the  diatoms,  i.e.,  ~1.8  per  day  at 
30  °C,  but  a  light  saturation  intensity  of  275  pmol  photons  m~2  s_1, 
reflecting  their  sun-adaptation  as  neutrally  buoyant  particles  at  the 
sea  surface.  In  addition,  both  of  these  functional  groups  utilized  only 
inorganic  forms  of  N  and  P.  Therefore,  selective  grazing  losses 
determined  the  outcome  of  the  competition  between  these  two 
fast-growing  groups  in  the  model,  unless  Si-limitation  prevailed. 


2.2.3.  Herbivory 

When  diatoms  were  present,  copepods  actively  selected 
this  food  item  in  the  model,  as  observed  on  the  WFS  (Kleppel 
et  al.,  1996).  instead,  the  mean  grazing  loss  of  the  smaller 
microflagellates  to  protozoan  micro-zooplankton  in  the  COM 
was  ~82%  of  their  realized  growth  rate  (Fahnenstiel  et  al., 
1995),  while  the  larger  diatoms  suffered  minimal  predation  by 
ciliates  (Strom  and  Strom,  1996).  Therefore,  we  imposed 
grazing  stresses  on  first  the  diatoms  and  microflagellates, 
with  minimal  losses  to  K.  brevis  (Turner  and  Tester,  1997) 
and  Trichodesmium  (O’Neil  and  Roman,  1992;  Guo  and  Tester, 
1994).  If  the  daily  food  demand  of  the  copepods  was  not  met 
by  diatoms,  they  then  consumed  protozoans  in  the  model 
(Walsh  et  al.,  2001b).  In  turn,  microflagellate  losses  to  proto¬ 
zoans  were  a  function  of  just  the  algal  prey  biomass  (Table  1), 
since  the  unicellular  herbivores  were  an  implicit  state  variable 
of  the  model. 

Accordingly,  the  model’s  total  phytoplankton  community  was 
grazed  by  both  copepods  and  protozoans,  which  were  eaten  in 
turn  by  competing  chaetognath  rivals  of  the  larval  clupeids 
(Fig.  2).  These  primary  carnivores  thus  consumed  herbivorous 
zooplankton  and  were  the  prey  of  higher  trophic  levels  (Fig.  2). 
Although  copepods  generally  shun  K.  brevis  (Tester  et  al.,  2000), 
the  dominant  smaller  WFS  calanoid  copepods,  e.g.,  Paracalanus 
quasimodo,  would  eat  this  species  if  other  prey  were  not  available 
within  downstream  North  Carolina  shelf  waters  (Turner  and 
Tester,  1997). 

We  also  found  cells  of  K.  brevis  in  the  mouths  of  femora 
turbinada  and  Centropages  velificatus  on  the  WFS  (Walsh  et  al., 
2003),  such  that  the  model’s  grazing  coefficient  was  0.01  on  K. 
brevis  (Table  1 ),  compared  to  0.70  on  diatoms  (Walsh  et  al.,  2003). 
However,  the  copepods  ate  mainly  diatoms  (Kleppel  et  al.,  1996) 
as  they  were  the  dominant  phytoplankton  community  on  the 
WFS  in  early  spring  of  the  “real  world”  (Gilbes  et  al.,  1996,  2002) 
and  the  model. 


2.2.4.  Autotrophic  processes 

In  contrast,  K.  brevis  had  a  maximal  growth  rate  of  only 
~0.8  day-1  (Walsh  et  al.,  2001a)  at  30  °C  in  our  model,  with  a 
larger  cellular  chlorophyll  content  of  ~5.0  x  10~6  pg  chi  cell-1 
and  a  light  saturation  intensity  of  65  pmol  photons  m~2  s_1 
(Table  1),  reflecting  their  shade-adaptation  (Walsh  et  al.,  2002). 
K.  brevis  can  utilize  urea,  vitamins,  amino  acids,  and  other  forms 
of  organic  nitrogen  at  varying  rates  (Baden  and  Mende,  1979; 
Shimizu  and  Wrensford,  1993),  providing  a  possible  niche  for 
these  and  other  dinoflagellates  (Dennison  and  Abal,  1999;  Glibert 
and  Terlizzi,  1999).  Although  the  half-saturation  constant  for  NH4 
uptake  by  K.  brevis  was  lower  than  that  for  urea  and  nitrate 
(Table  1),  uptake  of  all  three  nitrogen  pools  (labile  DON,  NH4,  and 
N03)  occurred  in  the  model. 


The  toxic  dinoflagellates  swam  vertically  at  a  speed  of 
1  mlr1,  in  response  to  a  light  cue  (Walsh  et  al.,  2002;  Kerfoot 
et  al.,  2004).  The  diatoms  aggregated,  settling  as  a  function  of 
their  biomass  (Walsh  and  Dieterle,  1994),  and  the  microflagellates 
did  not  sink.  The  diazotrophs  also  vertically  migrate  in  relation  to 
the  light  field  by  controlling  their  buoyancy  (Walsby,  1978,  1992; 
Villareal  and  Carpenter,  1990).  Thus,  vertical  ascents  and  descents 
of  Trichodesmium  in  the  model  occurred  at  a  speed  of  ~3.2  m  h_1 
(Lenes  et  al.,  2008). 

Finally,  Trichodesmium  had  a  similarly  slow  maximal  growth 
rate  of  ~0.7  day-1  at  30  °C  as  K.  brevis  (Lenes  et  al.,  2005,  2008). 
However,  the  nitrogen-fixers  were  assigned  a  light  saturation 
intensity  (/sat)  of  400  pmol  photons  m~2  s_1  (Rueter  et  al.,  1990, 
1992)  in  the  model,  reflecting  their  light-adaptation  (Lenes  et  al., 
2005).  Since  the  ambient  dissolved  nitrogen  gas,  N2,  does  not 
limit  growth  of  these  diazotrophs  on  the  WFS,  the  realized  gross 
growth  rate  for  Trichodesmium,  p3,  was  determined  by  the  least 
available  resource  of  light  or  nutrients  (Lenes  et  al„  2008),  where 
dFe  was  dissolved  iron: 


/t3  =  p'  mm 


P04 


OOP  \ 


dFe 


IsatO1  h/'sat'  \(!<po4+P04)’  (koop  +  DOP))  ’  (krfpe  +  dFe) 


(1) 


In  eq.  (1),  pT  was  the  maximum  gross  growth  rate,  adjusted  for 
temperature-stress,  T,  from  the  measured  temperature  profiles, 
assuming  that  the  Trichodesmium  growth  rate  was  maximal  at 
temperatures  >  20  °C  (Carpenter,  1983).  The  Michaelis-Menten 
half-saturation  constants  for  uptake  of  phosphate  (kP04),  dis¬ 
solved  organic  phosphorus  (kDop).  and  iron  (kdFe)  were  specific 
to  Trichodesmium.  Phosphorus  availability  was  calculated  as  TDP, 
where  preferential  uptake  was  given  to  P04  over  DOP  by  imposi¬ 
tion  of  a  smaller  half-saturation  constant  (Table  1). 

High  levels  of  alkaline  phosphatase  activity  have  been  mea¬ 
sured  in  phosphate-depleted  diazotroph  colonies  (Yentsch  et  al., 
1972;  Stihl  et  al.,  2001),  indicating  that  Trichodesmium,  like 
Karenia,  can  assimilate  DOP  (Vargo  and  Howard-Shamblott, 
1990).  Each  PFT  was  thus  assigned  a  unique  mean  C:N:Si:P:Fe 
molar  ratio  to  accentuate  their  competition  (Table  1).  In  conjunc¬ 
tion  with  different  half-saturation  constants,  the  ability  of  each 
PFT  to  utilize  different  nutrient  pools,  and  element-specific 
remineralization  rates,  helped  compensate  for  the  mathematical 
restrictions  associated  with  the  selection  of  the  Monod  limitation 
equations  and  fixed  species-specific  nutrient  ratios. 

We  computed  the  2001  daily  influx  of  iron  to  WFS  surface 
waters  as  perturbation  conditions  (Fig.  4).  Measured  atmospheric 
dust  concentrations  by  Joe  Prospero  at  Miami  and  daily  rainfall  at 
Tampa  International  Airport  were  used  to  calculate  both  wet  and 
dry  deposition  (Lenes  et  al.,  2008).  The  cellular  N/Fe  molar  ratio  of 
Trichodesmium  (Table  1)  was  1500  (Sanudo-Wilhelmy  et  al., 
2001).  Diatoms  required  10-fold  less  dFe  during  synthesis  of 
nitrate  reductase  for  nitrate  assimilation  in  the  model,  compared 
to  what  the  diazotrophs  needed  in  formation  of  nitrogenase 
during  nitrogen  fixation.  Therefore,  dFe  inputs  (Fig.  4)  mainly 
influenced  the  growth  of  diazotrophs,  via  Eq.  (1 ),  due  to  their  high 
Fe  demand. 

Trichodesmium  populations  of  the  model  were  minimally 
grazed,  but  they  did  undergo  cell  lysis  as  a  function  of  nutrient 
limitation  (Lenes  et  al.,  2005,  2008),  to  mimic  programmed  cell 
death  in  response  to  nutrient  stress  (Berman-Frank  et  al.,  2004, 
2007).  Photo-inhibition  of  Trichodesmium  was  ignored  in  HABSIM, 
because  high  saturation  light  intensities  (Carpenter  and 
Roenneberg,  1995)  allow  them  to  thrive  within  surface  summer 
waters  of  the  WFS  (Lenes  and  Heil,  2010),  while  Subramanian 
et  al.  (1999)  has  shown  that  Trichodesmium  both  releases  excess 
light  energy  as  fluorescence  and  possesses  high  levels  of  photo- 
protective  pigments. 
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2.2.5.  Bacteria 

Bacterial  secondary  production  can  be  an  order  of  magnitude 
higher  within  a  red  tide  on  the  WFS  than  in  surrounding  waters 
(Heil  et  al.,  2004).  Thus,  heterotrophs  may  be  competing  against 
K.  brevis  for  utilization  of  the  DON  (Mulholland  et  al.,  2004,  2006). 
Since  bacteria  also  compete  against  the  phytoplankton  for  both 
phosphate  and  ammonium  in  the  GOM  (Pakulski  et  al.,  2000), 
they  were  included  as  explicit  state  variables  (Table  1 ).  The  two 
distinct  functional  groups  of  the  model  were  ammonifying  and 
nitrifying  bacteria  (Darrow  et  al.,  2003;  Lenes  et  al.,  2008).  Like 
the  autotrophic  rates  of  the  phytoplankton,  the  heterotrophic 
processes  of  the  bacteria  were  impacted  by  the  temperature  fields 
(Kirchman,  2000;  Ducklow,  2000).  Their  temperature-regulated 
net  growth  rate  was  a  minimal  function  of  the  available  nitrogen 
or  phosphorus,  in  both  organic  and  inorganic  forms  (del  Giorgio 
and  Cole,  2000).  Specifically,  ammonifiers  assimilated  the  organic 
nutrient  pools  while  the  nitrifiers  utilized  the  inorganic  ones,  but 
in  N/P  ratios  of  10  in  the  model,  not  4  of  the  “real  world”.  Bacteria 
had  no  simulated  settling  velocities  in  the  model,  because  of  their 
small  size. 

On  the  WFS,  dissolved  organic  matter  (DOM)  tends  towards 
low  phosphorus  content,  with  a  large  molar  C/P  ratio  of  250 
(Walsh  et  al.,  2006).  Instead,  bacteria  have  a  smaller  molar  C/P 
ratio  of  53.  The  net  result  of  these  calculations  in  our  model  was 
that  the  bacteria  competed  with  the  phytoplankton  for  phos¬ 
phate,  when  labile  DOC  was  abundant.  This  result  was  consistent 
with  the  observations  that  DOC-rich  aquatic  systems  tend 
towards  P-limitation  (Kirchman,  1994).  Inclusion  of  two  bacterial 
pools  allowed  for  a  differential  response  in  nutrient  remineraliza¬ 
tion  dependent  upon  the  nutrient  concentrations  within  available 
pools. 

2.2.6.  Dead  fish 

Adult  and  juvenile  fishes  range  over  the  WFS,  harvesting  local 
zooplankton  food  supplies,  derived  in  turn  from  copepod  herbi¬ 
vores  eating  diatoms  and  from  copepod  omnivores  eating  proto¬ 
zoans  after  their  ingestion  of  microflagellates  (Houde  and  Chitty, 
1976;  Kleppel  et  al.,  1996).  Diazotrophs  (Guo  and  Tester,  1994) 
and  K.  brevis  (Turner  and  Tester,  1997)  are  instead  not  readily 
consumed  by  zooplankton,  nor  indirectly  by  teleosts.  Indeed,  K. 
brevis  may  harvest  the  fish  within  an  ecological  role  reversal  of 
direct  brevetoxic  poisonings  (Walsh  et  al.,  2006;  2009).  As  with 
surface  aggregations  of  diazotrophs,  positively  buoyant  dead  fish 
bloat  and  accumulate  at  frontal  boundaries  of  the  sea  surface 
(Gunter  et  al„  1948).  Such  dead  fish  are  thus  a  concentrated 
source  of  recycled  nutrients  within  the  euphotic  zone,  but  derived 
from  a  larger  area  than  that  of  the  local  water  column. 

These  fish  supplies  of  nutrients  may  be  equivalent  to  those 
derived  from  dying  Trichodesmium  at  the  beginning  of  a  red  tide 
within  offshore  waters  (Lenes  and  Heil,  2010).  Thus,  we  have  now 
included  within  the  HABSIM  model  (Fig.  2)  a  positive  feedback  of 
the  recycled  organic  nutrients  (DON  and  DOP)  from  decaying  fish, 
killed  by  K.  brevis.  Note  that  dissolved  silicon  was  not  replenished 
during  this  simulated  remineralization  process  of  nitrogen  and 
phosphorus. 

Yet,  zooplankton  biomasses  of  herbivores  and  omnivores  were 
not  explicit  state  variables  of  the  model.  Ingestion  demands  of 
two  fish  life  history  stages,  or  size  classes,  were  instead  used  to 
compute  the  cumulative  equivalent  nutrient  stocks  of  their 
zooplankton  prey,  ZP,  i.e.,  ZP^  2  of  Eq.  (2)  and  Table  1. 
d^  =  -8^+060[a{1_Ed)(ydpd) 

-Hl-EfcXrA)]  [1-0.30-0-0.01]  (2) 

The  first  pool,  ZP^,  represented  food  eaten  by  younger 
zoophagous  fishes,  while  raised  from  copepods,  grown  in  turn 


on  diatom  and  protozoan  prey.  We  assumed  that  most  of  these 
zooplanktivore  juvenile  and  larval  fish  would  be  young-of-the- 
year  clupeids,  undergoing  larval  recruitment  of  the  dominant 
spring  spawners  (Houde  et  al.,  1979),  i.e.,  the  thread  herring 
(O.  oglinum')  and  the  scaled  sardine  (Harengula  jaguana). 

The  second  group  of  fish  food,  ZP2,  represented  a  background 
pool  of  cumulative  prey  eaten  by  older  spawning  adult  fishes,  and 
released  when  these  fishes  were  poisoned.  The  adult  fish  pool  was 
set  as  a  homogeneous  nutrient  field  (Table  1),  provided  from 
mortality  of  dead  fish  during  toxic  events.  Therefore,  ZP2  was  set 
to  5.0  pmol  PN  kg-1  of  fish  body  nitrogen,  sampled  during  Octo¬ 
ber  2001  on  the  WFS  as  the  equivalent  amount  of  recycled 
ammonium  stocks  found  adjacent  to  dead  fish. 

To  specify,  when  and  where  in  the  model,  these  life  history 
stages  of  fishes  were  poisoned  by  toxic  levels  of  K.  brevis,  an  IF 
statement  was  employed.  When  the  concurrent  accumulation  of 
K.  brevis  reached  toxic  levels  (Landsberg,  2002),  the  cumulative 
supply  of  the  zooplankton  prey  assimilated  by  clupeoids, 
ZP=ZPi+ZP2,  were  released  as  DOM  at  that  grid  point.  The  rate 
of  nutrient  release  was  set  at  a  maximal  50%  of  ZP  per  day, 
reflecting  the  highest  dissolution  rate  of  dead  fish  in  warm  Florida 
waters  (Stevenson  and  Childers,  2004). 

Sloppy  feeding  by  copepods  of  this  simple  food  chain  leading 
to  younger  clupeiods  was  set  to  50%  (Walsh  et  al.,  2003),  such  that 
half  of  the  grazed  phytoplankton  became  a  DOM  source  within  eq. 
(2)  of  the  model.  Accordingly,  the  cumulative  ZP!  food  required 
by  early  life  history  stages  of  the  clupeoids  over  the  study  site 
(Fig.  1)  was  estimated  from  other  terms  of  Eq.  (2).  Of  these,  the 
first  one  of  the  settling  velocity  (wz)  was  set  to  zero,  in  relation  to 
the  vertical  coordinate  system  (er).  The  smaller  dead  fish  were 
also  assumed  to  float  at  the  surface,  like  the  adults,  as  a 
consequence  of  internal  gases  released  during  decomposition. 

In  the  second  term  of  Eq.  (2),  0.60  was  the  assimilation 
efficiency  of  copepod  herbivores,  feeding  poorly  on  just  diatoms, 
defined  by  the  sloppy  grazing  coefficient,  a,  the  specific  grazing 
rate,  yd,  and  a  respiration  rate  of  ed  (Walsh  et  al.,  2003).  When  the 
copepods  ate  protozoans,  the  efficiency  of  0.60  was  replaced  by 
one  of  0.80,  reflecting  higher  quality  food.  The  diatoms  were  a 
more  palatable  food  source  than  K.  brevis,  which  were  ingested  at 
a  much  smaller  rate,  yk,  with  a  respiration  rate  of  ek  (Walsh  et  al., 
2003).  The  daily  loss  components  of  the  fish  size  class  were 
calculated  in  the  last  term  of  Eq.  (2).  Respiration  loss  was  30%, 
while  reproduction  costs  to  the  larval  and  juvenile  fish  were  0%, 
and  removal  by  fish  piscivores  (groupers,  mackerels,  snappers) 
and  humans  was  1%. 


2.2.7.  Benthos 

On  the  sea  floor,  any  toxic  dinoflagellates  would  have  compe¬ 
titors  in  the  form  of  both  near-bottom  diatoms  and  benthic 
microflora  (Fig.  2).  From  prior  measurements,  the  chlorophyll 
biomass  of  photosynthetic  organisms  within  the  upper  one 
centimeter  of  WFS  sediments  was  greater  than  that  of  the  over- 
lying  water  column.  These  phytobenthos  were  likely  to  intercept 
any  ground  water  supplies  of  nutrients  (Hu  et  al.,  2006).  However, 
upon  recycling  of  their  body  N  and  P  by  bacteria,  like  pelagic 
diatom  sources  (Grill  and  Richards,  1964;  Darrow  et  al.,  2003), 
they  represent  another  source  of  low  N/P  nutrients  for  utilization 
by  the  diazotrophs  and  red  tides  on  the  WFS.  Accordingly,  the 
sediment  microflora  were  an  additional  state  variable  of  the 
model  (Fig.  2)  in  all  four  cases. 

We  assumed  that  the  benthic  microflora  were  diatoms 
(Darrow  et  al.,  2003),  with  a  maximum  growth  rate  of  ~0.5  day-1 
at  near-bottom  temperatures  of  10  DC.  Like  K.  brevis,  they  were 
also  shade-adapted,  with  again  a  light  saturation  intensity  of 
65  pmol  photons  m~2  s_1  (Table  1;  Darrow,  2008).  However,  unlike 
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Trichodesmium  and  the  heterotrophic  bacteria  of  the  model,  the 
benthic  microflora  only  used  inorganic  forms  of  the  nutrients  (P04, 
NH4,  N03,  and  Si04)  with  appropriate  values  of  their  half-saturation 
constants  drawn  from  the  literature  (Darrow  et  al.,  2003).  Based 
upon  temporal  changes  of  microfloral  biomass  between  April  and 
July  on  the  WFS,  as  well  as  between  April  and  June  on  the  Georgia 
shelf  (Nelson  et  al.,  1999),  a  grazing  loss  to  the  meiobenthos  was 
imposed  on  the  benthic  diatoms  of  the  model,  similar  to  that  of  the 
pelagic  diatoms  from  copepods.  In  the  thin,  well  oxygenated 
sediment  layer  of  the  model,  we  ignored  denitrification.  The  two 
benthic  bacterial  functional  groups  were  identical  to  those  in  the 
water  column. 

2.3.  Simulation  cases 

Four  cases  of  the  model  were  run  to  test  the  relative  role 
of  each  major  nutrient  source  during  2001.  The  baseline  of 
Case  1  included  three  nutrient  sources:  discharge  from  the  northern 
rivers:  Trichodesmium ;  and  dead  fish.  In  each  of  the  three  subse¬ 
quent  cases,  one  source  was  deleted  to  quantify  the  impact  on  the 
simulated  population  levels  of  K.  brevis,  compared  to  field  observa¬ 
tions  during  2001.  Case  2  had  no  river  pulse,  while  Case  3  had  no 
Trichodesmium  populations,  and  Case  4  had  no  fish  kills.  In  terms  of 
statistical  goodness-of-fits  of  the  model  cases  with  observed  cell 
counts  of  K.  brevis,  the  simulated  and  observed  chronologies  were 
compared  over  two  time  periods:  1  April-31  December  2001  and  1 
April-5  October  2001.  The  latter  results  are  designated  as  Cases  la, 
2a,  3a,  and  4a. 

2.4.  Validation  data 

In  this  study,  we  utilized  monthly  data  collected  during  three 
cross-shelf  sections  off  Tampa  Bay,  Sarasota,  and  Charlotte  Harbor 
and  two  along  shore  transects  from  1998-2002  (Fig.  1)  as  part  of 
the  NOAA/EPA  ECOHAB:  Florida  and  ONR  HyCODE  studies  (Walsh 
and  Kirkpatrick,  2008),  with  emphasis  here  on  the  2001  observa¬ 
tions.  Hydrographic  data  and  water  samples  were  collected  at  63 
locations,  with  a  Seabird  CTD  mounted  on  a  rosette  sampler  of 
Niskin  bottles.  Species-specific  phytoplankton  cell  counts  were 
performed  through  microscopic  examination,  both  during  the 
cruises  and  along  the  shoreline.  Vargo  et  al.  (2008)  have  described 
the  methodology  used  to  determine  some  of  the  inorganic  (N03, 
P04,  Si04)  and  total  dissolved  nutrients  (TDN,  TDP),  as  well  as  the 
chlorophyll  stocks. 

Winter  values  were  estimated  from  these  data  sources  and 
linearly  interpolated  to  initialize  the  model  (Table  2).  The  seaso¬ 
nal  mean  distributions  of  the  inorganic  nutrients  were  compiled, 
when  available  at  the  10-m  and  25-m  isobaths.  However,  con¬ 
tamination  of  the  2001  ammonium  samples  prevented  back 


calculation  of  DON  by  subtracting  nitrate,  nitrite  and  ammonium 
from  TDN.  Therefore,  DON  observations  were  not  available  for 
comparison  with  model  hindcasts  during  the  2001  simulations. 
We  instead  used  measurements  of  TDN  (TDN  =  N03  +  NH4+D0N) 
and  TDP  (TDP=P04+D0P)  for  nutrient  validation  of  the  model’s 
fidelity.  Finally,  sensitivity  simulations  were  run  both  to  derive 
initial  values  for  state  variables  that  lacked  in  situ  measurements 
and  to  evaluate  the  validity  of  the  model’s  solutions. 

In  terms  of  biomass  predictions  of  the  different  cases  of  the 
model,  validation  data  of  K.  brevis  cell  counts  during  2001  were 
obtained  from  the  Florida  HAB  historical  database  maintained  by 
the  Florida  Wildlife  Commission’s  Fish  and  Wildlife  Research 
Institute  (FWC  FWRI).  Coastal  circulation  dramatically  affects 
horizontal  position  of  K.  brevis  blooms  on  the  WFS  (Weisberg 
et  al.,  2009;  Walsh  et  al„  2009).  In  addition,  most  of  the  K.  brevis 
data  was  collected  in  “  an  event  response  mode”,  precluding  a 
consistent  temporal  time  series  at  any  one  location. 

However,  maximal  concentrations  of  phytoplankton,  detected 
either  at  the  sea  surface  by  satellites  (Walsh  et  al.,  2011a),  or  at 
depth  by  in  situ  sampling  of  cell  counts,  all  reflect  a  net  “patchy” 
realization  of  the  local  population  balances  of  the  growth  pro¬ 
cesses  of  photosynthesis  and  nutrient  uptake,  in  relation  to  the 
loss  processes  of  catabolism,  sinking,  and  grazing.  If  one  of  these 
factors  is  altered  significantly,  i.e.,  either  addition  of  nutrients,  or 
reduction  of  grazing  pressure,  the  next  set  of  values  of  temporal- 
spatial  maxima  will  accordingly  reflect  the  new  balances  of 
population  gains  and  losses  within  a  local  region  of  the  WFS 
(Fig.  1). 

Accordingly,  our  1-day  vertical  model  at  a  fixed  position 
(Fig.  1 )  represented  a  cumulative  WFS  catchment  basin  of  onshore 
and  alongshore  advection  of  phytoplankton  populations  during 
net  biological  growth  and  loss  processes  (Walsh  et  al.,  2003). 
Thus,  the  weekly  maximum  K.  brevis  cellular  concentrations  were 
extracted  from  the  FWC  FWRI  database,  over  both  the  surface 
(  <  3  m)  and  bottom  (within  5  m  of  bottom)  within  35  km  of  the 
West  Florida  coast  between  Cedar  Key  and  Naples  (Fig.  1),  for 
evaluation  of  the  sensitivity  cases  of  the  model. 

Such  an  approach  allowed  us  to  compare  the  simulated 
biochemical  bloom  dynamics  independent  of  the  horizontal 
physical  processes  and  sampling  bias,  since  this  region  (Fig.  1) 
comprised  >  80%  of  the  K.  brevis  samples  (Young  and  Christman, 
2006;  Heil  and  Steidinger,  2008).  These  maximal  K.  brevis  cell 
counts  were  converted  to  mmolCm~3  for  comparison  to  the 
model  output,  using  5.0  x  10~6  pgchl  cell  1  and  a  C:chl  ratio  of 
30  (Table  1).  Given  the  prior  concurrence  of  weekly  mean  cell 
counts  of  the  K.  brevis,  when  present  and  dominant  along  the 
West  Florida  coast,  with  satellite  estimates  of  the  total  phyto¬ 
plankton  community  (Walsh  et  al.,  2011a),  we  feel  that  these 
additional  observations  of  the  time  history  of  the  red  tides  during 


Table  2 

Validation  statistics  and  metrics  applied  to  the  Karenia  brevis  functional  group  by  case.  Statistics  for  cases  1  -4  are  calculated  based  on  the  full  temporal  model  simulations. 
Statistics  for  cases  la,  2a,  3a  and  4a  are  through  5  October.  Bold  values  indicate  the  “best  fit". 


Statistic 

Case  1 

Case  la 

Case  2 

Case  2a 

Case  3 

Case  3a 

Case  4 

Case  4a 

AE 

44.254 

5.475 

34.760 

-19.644 

-37.838 

-57.490 

-36.043 

-54.150 

AAE 

57.949 

24.120 

61.635 

27.712 

37.880 

57.586 

37.297 

56.575 

RMSE 

92.525 

42.655 

98.329 

54.252 

84.915 

99.904 

83.872 

98.091 

RMSE  IP 

2.379 

0.729 

2.528 

0.927 

2.183 

1.707 

2.156 

1.676 

MEF 

-0.462 

0.747 

-0.651 

0.590 

-0.231 

-0.390 

-0.201 

-0.340 

r2 

0.281 

0.745 

0.140 

0.631 

0.274 

0.634 

0.194 

0.151 

Average  error  (bias):  AE=(E(Pi-0 ,))/n;  average  absolute  error:  AAf— (Z!P,-  0,|)/n;  root  mean  squared  error:  RMSE=((S(P,  —  0,)2)/n)°'5:  general  standard  deviation:  RMSE/ 
P,  where  P=(ZP,)/n;  modeling  efficiency:  MEF=(X(0(  O)2  X(P;  0,)^)/X( Oj  O)2;  coefficient  of  determination:  r2= [SJOj — O)  ■  (P, — P)/(S(0,  —  O)2  •  S(P, — P)2)0'5]2,  where 
n=the  number  of  observations,  Oi=the  ith  of  n  observations,  P,=the  ith  of  n  predictions,  O  and  P  are  the  observation  and  prediction  averages  (statistical  equations  after 
Biber  et  al.,  2004;  Stow  et  al„  2009). 
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July-December  2001  allow  ease  of  interpretation  of  the  results  of 
the  four  case  studies,  from  this  admittedly  complex  model  of 
Florida  red  tides  of  K.  brevis. 

2.5.  Sensitivity  analyses 

Testing  the  reliability  of  an  ecological  model  requires  assess¬ 
ment  of  the  input  parameters.  A  sensitivity  analysis  of  these 
ecological  parameters  also  provides  description  of  the  needs  for 
model  improvement,  indicating  which  state  variables  have  the 
greatest  impact  on  the  ability  of  the  model  to  replicate  observa¬ 
tions.  In  this  model,  we  tested  the  sensitivity  of  95  parameters  in 
relation  to  the  computed  values  of  K.  brevis,  in  statistical  compar¬ 
isons  with  the  observations.  The  sensitivity  of  each  model  para¬ 
meter  was  calculated  (Franks  et  al.,  1986;  Fasham  et  al.,  1990; 
Ji  et  al.,  2006)  by  Eq.  (3),  in  which  A  Parameter  was  the  amount  of 
variation  of  a  1%  departure  from  the  baseline  value. 

S=  _ Af/F _  (3) 

A  Parameter/Parameter 

In  Eq.  (3),  S  was  the  index  of  sensitivity  and  AF  was  the  change 
in  concentration  of  K.  brevis,  compared  to  its  value  in  the  baseline 
simulation,  F,  caused  by  fluctuation  of  the  specified  model 
parameter.  The  prior  studies  defined  a  parameter  to  be  sensitive, 
if  S  >  0.5.  The  sensitivities  of  the  simulated  K.  brevis  blooms  to 
some  of  the  model  parameters  are  shown  in  Fig.  5.  Only  the 
parameters  with  an  index  of  >  0.1  are  considered  here  (19  out  of 
95  parameters).  We  shall  find  that  of  the  model’s  95  parameters, 
only  three  were  considered  significantly  sensitive:  K.  brevis 
grazing  (y4);  the  molar  phosphorus  to  carbon  ratio  of  K.  brevis 
((C:P)4);  and  the  maximum  growth  rate  of  Trichodesmium  (p3). 

Six  different  statistical  criteria  (Table  2)  were  evaluated  for 
each  of  the  four  separate  cases  to  assess  the  predictive  ability  and 
efficiency  of  this  version  of  the  HABSIM  model.  The  model  results 
were  compared  against  weekly  surface  and  near-bottom  maximal 
cell  counts,  when  available,  yielding  estimates  of:  average  error 
(AE);  average  absolute  error  (AAE);  root  mean  squared  error 
(RMSE);  general  standard  deviation  (RMSE/P);  modeling  effi¬ 
ciency  (MEF);  and  the  coefficient  of  determination,  i.e.,  the  square 
of  the  correlation  coefficient  (r2).  We  utilized  this  set  of  quanti¬ 
tative  metrics  for  univariate  comparison  of  model  results  with 
observations  as  outlined  by  Biber  et  al.  (2004)  and  Stow  et  al. 
(2009):  average  error  (AE),  average  absolute  error  (AAE),  root 
mean  squared  error  (RMSE),  general  standard  deviation  (RMSE/P), 
modeling  efficiency  (MEF)  and  correlation  coefficient  (r2).  The 
coefficient  of  determination  measures  the  tendency  of  the  pre¬ 
dicted  and  observed  values  to  vary  together,  with  an  ideal  value  of 
1.  The  root  mean  square  error,  average  error,  average  absolute 
error,  and  general  standard  deviation  all  calculate  the  size  of  the 
difference  between  predicted  and  observed  values,  where  values 
near  zero  indicate  a  close  match.  The  modeling  efficiency  repre¬ 
sents  how  well  a  model  predicts  relative  to  the  average  of  the 
observations,  where  values  near  one  indicate  a  close  match 
between  observations  and  predictions  and  a  value  less  than  zero 
indicates  the  observation  average  is  a  better  predictor  than  the 
model  results.  The  equations  for  each  metric  were  listed  in 
Table  2  and  a  more  comprehensive  description  of  the  use  of  these 
metrics  is  available  in  Stow  et  al.  (2009). 


3.  Results 

3.1.  Phytoplankton 

The  phytoplankton  community  of  the  model’s  baseline  case 
reproduced  initiation  of  the  spring  transition  of  microalgae,  with 
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Fig.  5.  Sensitivity  of  the  95  model  parameters,  of  which  only  those  parameters 
with  a  sensitivity  of  >  0.1  are  shown. 


a  diatom  pulse  by  10  April  of  >  1  pmolCkg-1  throughout  the 
water  column  (Fig.  6(a)).  These  fast  growing  siliceous  phyto¬ 
plankton  outcompeted  the  other  phytoplankton  for  the  inorganic 
nutrients  injected  at  the  surface  on  16  March.  The  simulated 
injection  of  nitrate  represented  the  annual  pulse  of  nutrients  from 
the  northeastern  COM  estuaries  (Gilbes  et  al.,  1996,  2002),  with 
subsequent  surface  chlorophyll  plumes.  The  modeled  response 
was  a  ~2  week  delay,  relative  to  the  observed  chlorophyll  feature 
on  the  northern  WFS. 

Such  a  time  delay  was  consistent  with  an  estimated  transit 
time  from  the  nutrient  source  region  to  the  model  site.  The 
simulated  diatoms  were  quickly  removed  from  the  model’s  water 
column  by  grazing  and  sinking,  leading  to  a  near  bottom  accu¬ 
mulation  of  >  1  pmolCkg-1  during  April  (Fig.  6(a)).  Note  that 
the  prior  surface  phytoplankton  stocks  of  >  1  pg  chi  1~]  observed 
during  the  3-5  March  2001  ECOHAB  cruise  had  also  sunk  out  of 
the  water  column  near  the  model’s  location  by  the  3-6  April 
cruise  (Fig.  7). 

When  pelagic  diatoms  reached  the  benthic  layer  in  the  model, 
they  were  converted  to  detritus  and  subjected  to  benthic  reminer¬ 
alization.  This  nutrient  source  fed  an  increase  in  microflagellate 
production,  leading  to  a  maximum  biomass  of  >  1  pmol  C  kg-1  on 
18  April  in  the  model  (Fig.  5(b))  despite  high  grazing  rates  (Table  1). 

Upon  subsequent  atmospheric  additions  of  iron-rich  Saharan 
dust  during  June  (Fig.  4),  the  simulated  T.  erythraeum  populations 
increased  by  an  order  of  magnitude  within  the  water  column, 
reaching  a  subsurface  July-August  maximum  of  >40  pmol  C  kg-1 
or  ~  2.5  pg  chi  1  1  (Fig.  6(c)).  The  opportunistic,  fast  growing  diatom 
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population  then  assimilated  these  nutrients  released  by  the  diazo- 
trophs,  accumulating  a  biomass  of  >3pmolCkg_1  in  July 
(Fig.  6(a)).  Once  again,  as  observed  in  prior  seasonal  succession  of 
phytoplankton  on  the  WFS  (Walsh  et  al.,  2003),  microflagellates 
followed  in  the  model,  with  an  August  accumulation  of 
>  2  pmol  C  kg-1  (Fig.  6(b)),  but  only  ~0.4  pgchl  l-1. 

Nutrients  recycled  from  the  precursor  diazotroph  bloom 
(Fig.  6(c))  also  fueled  slow  growth  of  K.  brevis,  which  reached 
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Fig.  6.  The  simulated  daily  phytoplankton  concentration  (mmol  Cm3  or 
pmol  C  kg-1)  of  (A)  diatoms,  (B)  microflagellates,  (C)  Trichodesmium,  and  (D)  Karenia 
brevis  over  the  20-m  water  column.  Note  change  in  scales. 


specified  fish  killing  concentrations  of  >  4.5  pmol  C  kg-1,  or 
1 .8  pg  chi  1  \  by  late  July  (Fig.  6(d)).  The  simulated  K.  brevis 
population  initially  accumulated  biomass  near  the  bottom  of  the 
September  water  column  (Fig.  6(d)),  due  to  its  low  light  satura¬ 
tion  intensity  (Table  1 ).  This  bottom  development  persisted,  until 
sufficient  biomass  had  accumulated  to  self-shade  the  HAB  popu¬ 
lation.  At  this  point,  the  early  HAB  of  K.  brevis  migrated  to  the 
surface  layer  of  the  model,  both  concentrating  and  co-locating 
the  red  tide  with  positively  buoyant  dead  fish  by  early  October 
(Fig.  6(d)). 

Simulated  K.  brevis  stocks  reached  ~  12  pmol  C  kg  1  on  10 
August,  or  ~4.8pgchll~1  (Fig.  6(d)).  This  was  similar  to  the 
chlorophyll  biomass  of  the  red  tide  as  measured  at  this  time  on 
the  inner  WFS  during  2001  (Vargo  et  al.,  2004,  2008).  Earlier  on 
28  July  2001,  satellite  imagery  suggested  a  total  phytoplankton 
stock  of  8-10  pg  chi  1  1  at  mid-shelf,  compared  to  1-2  pg  chi  1  1 
on  3  July  2001  (Fig.  8). 

Accordingly,  this  satellite  time  sequence  was  consistent  with  the 
model’s  surface  estimates  of  the  whole  phytoplankton  community  of 
~7.4pgchlD1  on  10  August  2001,  composed  of  4.8  pg  chi  I  1  of 
K.  brevis,  1.4pgchll_1  of  Trichodesmium,  0.8pgchll_1  of  diatoms, 
and  0.4pgchll_1  of  microflagellates.  Subsequently,  K.  brevis  grew 
steadily  throughout  the  fall  on  nutrients  harvested  from  decaying  fish 
and  diazotrophs  of  the  model,  to  surface  stocks  of  >  200  pmol  C  kg-1 
in  September,  or  >  80  pg  chi  1  1  (Fig.  6(d)). 

Satellite  imagery  of  the  phytoplankton  community  on  3  July 
(Fig.  8(a))  and  28  July  2001  (Fig.  8(b))  provided  an  initial 
assessment  of  the  diazotrophs  and  dinoflagellates  at  the  mid / 
inner  shelf.  On  3  July,  Trichodesmium  dominated  the  phytoplank¬ 
ton  community.  By  28  July,  the  optical  signal  (Carder  et  al.,  2007; 
Cannizzaro  et  al.,  2008)  of  K.  brevis  had  become  co-dominant  with 
the  diazotrophs,  as  observed  during  earlier  red  tides  (Gunter  et  al., 
1948;  Walsh  and  Steidinger,  2001). 

Once  again,  this  successional  sequence  of  phytoplankton, 
inferred  from  space,  was  consistent  with  model  results.  In  response 
to  ambient  dFe  and  P  stocks  of  the  model  and  then  deposition  of 
Saharan  dust  during  June-July,  simulated  growth  of  Trichodesmium 
led  to  population  increments  of  ten-fold  above  background  levels  by 
late  June  2001  and  one  hundred  fold  by  late  July  (Fig.  6(c)),  as 
observed  then  on  the  WFS  (Lenes  and  Heil,  2010). 

Little  to  no  surface  populations  of  Trichodesmium  were  present 
within  shelf  waters  during  April  2001,  with  initial  detection  of 
elevated  biomass  on  the  WFS  during  15  June  2001.  A  maximum 
concentration  of  3.0  x  106  trichomes  1  1  was  found  on  15  July 
within  surface  waters  of  the  inner  shelf  (Lenes  and  Heil,  2010). 


A  B 


3-6  April  2001. 
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Fig.  8.  The  SeaWiFS  satellite  estimate  of  the  optical  properties  of  red  tides  and  diazotrophs  in  the  eastern  Gulf  of  Mexico  during  (A)  3  July  and  (B)  28  July  2001.  The 
reddish-black  color  in  panel  B  is  a  red  tide  of  Karenia  brevis  and  the  pale  blue  color  in  both  panels  is  co-occurring  Trichodesmium,  based  upon  a  backscatter  algorithm 
(Carder  et  al.,  2007;  Cannizzaro  et  al.,  2008).  (For  interpretation  of  the  references  to  color  in  this  figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


The  observed  maximum  population  was  partially  the  result  of 
frontal  aggregation  within  wind  rows.  We  were  thus  not  surprised 
at  the  model’s  underestimate  (Fig.  6(c))  from  the  1-D  simulations, 
given  the  present  HABSlM’s  inability  to  aggregate  biomass  horizon¬ 
tally  in  a  3-D  mode. 

Simulated  surface  K.  brevis  stocks  successfully  reproduced  the 
initiation  and  early  maintenance  of  the  observed  K.  brevis  blooms, 
or  weekly  maxima  (open  circles  of  Fig.  9(a)),  with  an  r2  of  0.75 
(Case  la  of  Table  2)  from  1  April-5  October  2001.  Comparison  of 
the  baseline  results  of  the  model  to  the  observed  weekly  maximal 
K.  brevis  concentrations  yielded  an  MEF  of  0.75  for  Case  la.  These 
results  demonstrated  the  model’s  ability  to  simulate  changes  of 
the  surface  biomass,  as  seen  by  satellites,  in  relation  to  vertical 
mixing  of  the  population  during  wind  events.  For  example,  the 
decline  in  the  surface  biomass  in  early  September  (Fig.  9(a))  was 
accompanied  by  an  increase  in  the  near  bottom  population 
(Fig.  9(b)),  not  a  decrease  in  the  total  integrated  biomass. 

Within  the  context  of  a  3-D  “real  world”,  however,  the  vertical 
dispersion  of  HABs  could  be  a  significant  factor  in  both  termina¬ 
tion  and/or  horizontal  spreading  of  the  bloom,  in  response  to 
variations  in  the  direction  of  transport  at  different  depths 
(Weisberg  and  He,  2003).  The  1-D  version  of  HABSIM  performed 
poorly  in  reproducing  the  termination  of  the  surface  Karenia  HABs 
after  15  October  2001  (Fig.  9(a)).  The  baseline  results  of  the  model 
over  the  whole  simulation  period,  from  April  to  December  2001, 
had  an  r2  of  only  0.28,  with  an  MEF  of  -0.46  (Case  1  of  Table  2). 
The  average  error  increased  from  5.5  of  Case  la  to  44.3  of  Case  1 
(Table  2). 

There  are  several  possibilities  to  explain  this  failure  of 
long-term  model  fidelity,  comprised  -  as  usual  -  of  unknown 
phytoplankton  loss  terms  (Walsh,  1983).  Previous  simulations 
have  indicated  that  advective  losses,  due  to  shifts  in  the 
prevailing  winds,  play  a  critical  role  in  transport  and  disper¬ 
sion  of  the  surface  K.  brevis  populations  away  from  the  coast 
(Weisberg  et  al.,  2009;  Walsh  et  al.,  2009).  Under  conditions  of 
thermal  stratification,  the  shelf-wide  upwelling/downwelling 
responses  to  synoptic  winds  are  asymmetrical  (Weisberg 
et  al.,  2001)  and  favor  the  shoreward  progression  and  net 
retention  of  material  properties  in  the  bottom  frictional 
boundary  layer  (Weisberg  et  al.,  2009).  It  is  also  apparent  that 
near-surface  material  tracers  persist  longer  coherent  patches 
or  satellite-observed  features  during  periods  of  thermal  stra¬ 
tification  (Jolliff  et  al.,  2003,  2008).  As  the  Gulf  of  Mexico 
basin-wide  thermal  stratification  regularly  breaks  down  in 
early  October  due  to  a  decline  in  solar  forcing  and  an  increase 
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Fig.  9.  The  simulated  daily  baseline  (Case  1)  Karenia  brevis  concentrations 
(mmol  C  m  :i ;  solid  line)  at  the  (A)  surface  and  (B)  near  bottom  layer  in  relation 
to  the  observed  weekly  maximum  Karenia  brevis  concentrations  (O)  within  35  km 
of  the  West  Florida  coast  between  Cedar  Key  and  Naples.  Observed  surface  cell 
counts  were  collected  at  depths  <  3  m  and  near  bottom  counts  were  collected 
between  15-20  m. 


in  episodic  latent  heat  flux  losses  (Jolliff  et  al.,  2008),  the 
retentive  asymmetry  effects  are  likely  diminished  and  more 
vigorous  vertical  mixing  ensues.  Thus  after  October  1,  physical 
dispersion  of  any  accumulated  biomass  patches  in  the  “real 
world”  were  not  captured  in  the  1-D  model. 

In  addition,  the  model  utilized  constant  loss  terms  of  1% 
grazing  and  1%  lysis  per  day  (Table  1).  It  has  been  suggested  that 
K.  brevis  is  capable  of  programmed  cell  death  under  nutrient 
stress,  viral  lysis,  and  attack  by  algicidal  bacteria  (Doucette  et  al., 
1999;  Paul  et  al.,  2002;  Roth  et  al.,  2008;  Van  Dolah  et  al.,  2009), 
although  the  rates  are  yet  to  be  quantified.  Therefore,  revisions 
of  the  loss  term  and  coupling  to  a  3-D  circulation  model 
could  compensate  for  the  model’s  inability  to  reproduce 
bloom  termination. 
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Dissolved  Nitrogen  (TDN=N03+NH4+D0N)  and  (B)  Total  Dissolved  Phosphorus 
(TDP— [J04  1  DOP)  in  relation  to  the  observed  TDN  and  TDP  at  station  29  (27.2'N, 
-82.9W)  as  measured  during  monthly  ECOHAB  cruises.  Open  circles  represent 
field  data  while  the  solid  lines  are  model  results. 


3.2.  Nutrients 

In  early  April  2001,  the  simulated  surface  TDN  concentration 
of  the  baseline  case  was  ~5.8  pmol  N  kg-1  (Fig.  10(a)),  somewhat 
elevated  from  earlier  estimates  of  the  mid-shelf  background  TDN 
concentrations  of  ~5  pmol  N  kg  1  (Walsh  et  al.,  2003).  This 
model  result  was  also  consistent  with  the  surface  TDN  concentra¬ 
tion  measured  at  station  29  on  the  3-6  April  2001  ECOHAB  cruise. 
We  attribute  these  initial  spring  stocks  of  dissolved  nitrogen  off 
central  Florida  to  nutrient  sources  from  mainly  the  northern  Gulf 
estuaries. 

Subsequently,  the  simulated  TDN  increased  2-fold  over  the 
summer  months  to  a  maximum  of  ~  11.2  pmol  N  kg-1  by  early 
August  2001  (Fig.  10(a)).  Such  an  increment  resulted  from  both  N2 
fixation  by  Trichodesmium  and  decay  of  poisoned  fish,  with  a  nitrogen 
isotopic  signature  attributed  to  trophic  transfers  from  nitrate-grown 
diatoms  at  base  of  their  separate  food  chain  (Walsh  et  al.,  2009). 
Accordingly,  the  recycled  DON  of  the  model’s  baseline  case  (Fig.  11(f)) 
was  responsible  for  the  majority  of  this  TDN  increase. 

The  DON  was  quickly  converted  to  NH4  (Fig.  11(b))  and  N03 
(Fig.  11(a))  through  bacterial  processes  of  the  model,  where  it 
became  bioavailable  for  all  components  of  the  phytoplankton 
community.  Increases  in  both  the  simulated  diatom  and  microfla- 
gellate  populations  (Fig.  6(a)  and  (b))  directly  resulted  from  this  N 
supply.  Thus,  both  the  observed  (open  symbols  of  Fig.  11)  and 
simulated  nitrate  stocks  were  usually  less  than  0.5  pmol  Nkg-1 
during  2001.  Although  ammonium  data  were  not  available  in  2001, 
the  simulated  seasonal  maximum  of  >  0.2  pmol  N  kg-1  during 
August  2001  (Fig.  11(b))  was  the  same  as  those  maximal  amounts 
observed  above  the  10-50  m  isobaths  of  the  WFS  during  August 
1998  and  1999  (Walsh  et  al.,  2003). 


Fig.  11.  The  simulated  daily  surface  concentration  ( — )  and  measured  monthly  isobath  averages  (111  =  10  m;  0=25  m)  of  (A)  N03,  (B)  NH4,  (C)  P04,  (D)  dFe,  (E)  Si04, 
(F)  DON,  (G)  DOP,  (H)  ammonifying  bacteria  (Bi),  and  (I)  nitrifying  bacteria  (B2).  Note  change  in  scales. 
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Given  the  lack  of  a  Si  sources  associated  with  the  decay 
of  Trichodesmium  and  dead  fish,  the  diatoms  quickly  depleted 
the  model’s  August  Si04  concentrations  to  <  0.5  pmol  Si  kg-1 
(Fig.  11(e)).  This  was  again  consistent  with  both  the  Si04  concentra¬ 
tions  measured  at  the  10-m  isobath,  and  to  a  lesser  effect,  at  the  25-m 
isobaths  during  August  2001  (Fig.  9(e)),  as  well  as  during  August  1998 
and  1999  (Walsh  et  al.,  2003).  Note  that  the  half-saturation  constant 
for  silica  uptake  by  the  diatoms  was  0.5  pmol  Si  kg-1  in  the  model 
(Table  1 ).  Therefore,  the  seasonal  Si-limitation  of  the  diatom  popula¬ 
tion,  together  with  the  high  grazing  rates  on  the  microflagellates  kept 
both  of  their  potentially  fast-growing  populations  in  check,  allowing 
K.  brevis  to  dominate  the  fall  phytoplankton  community  (Fig.  6), 
despite  the  10  x  slower  growth  rates  of  the  HABs. 

While  the  simulated  TDN  pattern  was  consistent  with  the 
measured  concentrations,  the  model  still  predicted  a  significant 
delay  of  phytoplankton  biomass  accumulations  (~5  weeks), 
relative  to  those  found  at  station  29.  This  discrepancy  resulted 
from  the  model’s  underestimate  of  spring-summer  stocks  of  DOP 
(Fig.  11(g)),  and  hence  the  slow  initial  growth  of  the  modeled 
Trichodesmium  population. 

Competition  for  available  P  among  the  other  phytoplankton 
and  bacterial  pools  of  the  model  had  reduced  the  simulated  DOP 
concentrations  by  25-50%,  relative  to  the  observed  concentrations  at 
the  10-m  and  25-m  isobaths  by  the  end  of  June  2001(Fig.  11(g)).  In 
contrast,  the  model  correctly  predicted  the  P04  at  the  end  of  that 
month  (Fig.  1 1(c)).  But,  since  DOP  dominated  the  TDP  pool  (Fig.  10), 
such  P-limitation  suppressed  Trichodesmium' s  already  slow  growth 


rate,  despite  sufficient  Fe  stocks  by  this  date  (Fig.  11(d)).  Of  course,  if 
the  model  had  instead  employed  a  N/P  regeneration  ratio  of  4, 
rather  than  16,  phosphorus  limitation  of  Trichodesmium  would  have 
been  much  less,  and  the  lag  time  of  biomass  accumulation  would 
have  only  been  ~  9  day.  Such  sensitivity  in  the  model  to  formulation 
of  the  “correct”  N/P  ratios  of  nutrient  regeneration  was  shown  in  the 
highest  sensitivity  index  (Fig.  5)  of  the  C/P  ratio  of  the  Karenia  HABs, 
followed  by  that  of  the  growth  rate,  fi3,  of  Trichodesmium  (Fig.  5). 


3.3.  Other  model  cases 

In  Case  2  of  no  river  nutrient  loadings,  the  model  still  predicted  a 
K.  brevis  HAB  with  a  surface  value  of  ~ 225  mmol  Cm~3  on  15 
September  2001  (Fig.  12(a)),  compared  to  ~325  mmol  C  on  the 
same  day  of  the  baseline  case  (Fig.  9(a)).  A  few  days  earlier,  a 
similar  amount  of  ~250  mmol  C  m  !  of  near-surface  dinoflagellate 
stock  had  been  observed  that  month  on  the  WFS  (Fig.  9(a)).  This 
reduction  of  only  30%  between  the  two  model  cases,  which 
bracketed  the  observations,  suggested  that  the  spring  nutrient 
inputs  from  the  northern  estuaries  and  remineralization  of  the 
subsequent  diatom  population  were  not  causative  agents  of  red 
tides  on  the  WFS.  These  model  results  instead  reflected  the 
secondaiy  role  of  dead  fish  produced  in  this  parallel  food  chain  for 
the  maintenance,  but  not  the  initiation  of  red  tides  (Walsh  et  al., 
2009).  This  decreased  model  fidelity  was  confirmed  by  an  overall  r2 
of  0.14  for  Case  2,  compared  to  one  of  0.28  for  Case  1  (Table  2). 


Fig.  12.  The  simulated  daily  Karenia  brevis  concentrations  (mmol  Cm3)  for  Case  2  (no  river  pulse)  within  the  (A)  surface  and  (B)  near  bottom  layers,  for  Case  3  (no 
Trichodesmium)  within  the  same  (C),  (D)  layers;  and  for  Case  4  (no  dead  fish)  within  the  same  (E),  (F)  layers.  Note  change  in  scales. 
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The  removal  of  N2  fixers  in  Case  3  led  to  an  insignificant  red  tide 
of  K.  brevis  on  15  October  2001  of  <  2.5  mmol  C  m  3  within  surface 
waters  of  the  model  (Fig.  12(c)).  Ambient  nutrient  stocks  supported 
slow  growth  of  K.  brevis,  reaching  a  maximum  surface  concentration 
of  ~4  mmol  Cm-3  by  early  November.  While  this  HAB  concentra¬ 
tion  would  have  been  large  enough  to  force  closures  of  commercial 
harvests  from  shellfish  beds,  the  red  tides  of  the  model  then  did  not 
exceed  those  levels  needed  to  cause  significant  fish  kills. 

Similarly,  the  removal  of  dead  fish  in  Case  4  restricted  the 
maximum  simulated  concentration  of  K.  brevis  to  ~2.7  mmol  C  m~3 
on  15  October  2001  (Fig.  12(e)).  After  this  date,  the  absence  of  a 
dead  fish  maintenance  supply  of  nutrients  led  to  an  accumulated 
HAB  of  only  0.8  mmol  C  m  3  by  31  December  2001  in  Case  3 
(Fig.  12(e)),  compared  to  1.5  mmol  C  m~3  on  the  same  day  in  case 
2  (Fig.  12(c)).  Recall  that  very  small  amounts  of  K.  brevis  were  also 
observed  near  the  sea  surface  on  that  day  (Fig.  9(a)),  but  presumably 
as  a  result  of  advective  and/or  grazing  losses,  not  in  response  to 
meager  supplies  of  recycled  nutrients.  The  last  sensitive  parameter 
of  the  model  was  that  of  grazing  stress,  A4,  on  K.  brevis  (Fig.  5),  which 
can  vary  tenfold  from  one  year  to  the  next  on  the  WFS  (Walsh  et  al., 
2011a). 


4.  Discussion 

In  a  shallow  WFS  water  column  of  20  m,  rapid  sinking  velocities 
of  diatom  blooms  (Walsh  and  Dieterle,  1994)  and  zooplankton  fecal 
pellets  (Walsh  et  al.,  2003)  ensured  arrival  of  phytodetritus  and 
egested  herbivore  debris  in  less  than  one  day.  Based  on  conventional 
wisdom,  the  grazing  fluxes,  AP,  varied  for  each  phytoplankton  state 
variable,  as  a  function  of  their  assessed  palatability  to  copepods  and 
protozoans.  When  herbivores  are  explicit  state  variables  of  a  model 
(Walsh,  1976),  they  provide  both  grazing  closure  and/or  linkage  to 
other  components  of  the  coastal  food  web. 

With  respect  to  grazing  on  K.  brevis  by  copepods,  the  crusta¬ 
ceans  are  thought  to  be  ineffective  herbivores  of  these  HABs  by 
observationalists  (Kleppel  et  al.,  1996;  Turner  and  Tester,  1997; 
Lester  et  al.,  2008).  Yet,  grazing  controls  are  exerted  by  calanoid 
copepod  herbivores  on  the  WFS  in  numerical  models,  albeit  on 
diatoms  not  toxic  dinoflagellates  (Milroy  et  al.,  2008).  Indeed,  the 
clupeid  fish,  which  K.  brevis  harvests  as  a  supplemental  nutrient 
source  (Walsh  et  al.,  2009),  are  supported  by  a  short  parallel  food 
chain,  based  in  turn  on  spring  blooms  of  diatoms,  mainly 
Skeletonmea  costatum  (Saunders  and  Glenn,  1969),  eaten  by 
calanoid  copepods.  Furthermore,  based  on  their  specific  rates  of 
grazing  (Dagg,  1995)  and  observed  abundances  during  September 
1998,  copepod  herbivores  may  then  have  consumed  as  much  as 
50%  of  the  daily  production  (Sutton  et  al.,  2001)  of  the  previous 
near-bottom  summer  diatom  stocks  of  3-4pgchll_1  on  the 
middle  WFS. 

By  late  fall,  however,  no  diatom  food  was  left  for  the 
surviving  copepods  (Walsh  et  al.,  2003).  These  facultative 
herbivores  must  instead  consume  first  protozoans  (Kleppel 
et  al.,  1996),  and  then  populations  of  K.  brevis  in  the  “real 
world”.  Indeed,  with  no  alternative  diatom  prey,  the  copepods 
will  eat  these  toxic  dinoflagellates  in  coastal  waters  (Turner 
and  Tester,  1997),  preventing  accumulation  of  ungrazed  toxic 
summer  red  tides  on  the  WFS  (Milroy  et  al.,  2008)  of 
>1  pgchll-1,  as  was  observed  in  1966,  1972,  1993,  1998, 
and  2000  (Walsh  et  al.,  2006). 

On  the  30-m  isobath  off  Sarasota  and  Tampa,  FL,  for 
example,  the  dominant  copepod  genera  were  Paracalanus, 
Temora,  Centropages,  Oithona,  Oncaea,  and  Corycaeus  in  both 
September  1998  (Sutton  et  al.,  2001)  and  1999  (Lester  et  al., 
2008).  Although  these  copepods  generally  shun  K.  brevis 
(Tester  et  al.,  2000),  the  dominant  smaller  WFS  calanoid 


copepods,  e.g.,  P.  quasimodo,  will  eat  this  HAB,  if  no 
other  prey  are  available,  within  downstream  North  Carolina 
shelf  waters  (Turner  and  Tester,  1997).  Moreover,  long-term 
changes  of  the  phytoplankton  abundances  and  species 
composition  on  the  WFS  have  been  attributed  to  reduced 
grazing  pressures,  after  increased  consumption  of  the  herbi¬ 
vores  by  zooplanktivores  during  trophic  cascades  (Walsh 
et  al.,  2011a). 

Yet,  these  parameters  of  algal  palatability  were  just  implicit 
fixed  rates  of  phytoplankton  abundances  in  the  present  1-D 
simulations.  No  consideration  was  given  to  the  explicit  popula¬ 
tion  dynamics  of  the  zooplankton  herbivores,  which  both  eat 
these  microalgae  and,  in  turn,  are  subject  to  predation  by  clupeid 
and  chaetognath  carnivores.  Proper  grazing  stress  formulations 
are  thus  required  in  future  models,  as  well  as  time-dependent 
nutrient  regeneration  formulations  and  algal  lysis,  to  accurately 
derive  species  succession  and  HAB  maintenance  within  the 
context  of  3-D  physical  dispersion  loss  terms  as  part  of  the 
mechanisms  driving  bloom  termination. 

The  Coastal  Ocean  Monitoring  System  (COMPS)  at  the  Uni¬ 
versity  of  South  Florida  (USF)  was  developed  as  part  of  the 
Regional  Coastal  Ocean  Observing  System  (RCOOS)  for  the 
southeastern  United  States  (Weisberg  et  al.,  2009).  The  past 
COMPS  consisted  of  tide  gauges,  surface  meteorological  sen¬ 
sors,  buoys  with  acoustic  Doppler  current  profilers  (ADCPs), 
temperature  (T)  and  salinity  (S)  sensors,  high  frequency  radar 
for  surface  current  mapping,  and  bottom  stationed  ocean 
profilers  (BSOPs)  for  discrete  measurements  of  T,  S,  and  chlor¬ 
ophyll  fluorescence  at  a  few  locations  (see  http://ocgweb.mar- 
ine.usf.edu).  Such  COMPS  data,  supplemented  with  glider 
observations  could  now  provide  both  the  offshore  boundary 
conditions  and  validation  criteria  for  real-time  improvements 
of  the  fidelity  of  future  operational  HAB  models.  With  revisions 
to  a  3-D  version  of  HABSIM,  including  time-dependent  N/P 
remineralization  ratios,  explicit  herbivore  state  variables,  and 
physical  dispersion  losses,  we  will  next  explore  the  reasons  for 
the  large  HAB  of  2001,  compared  to  the  small  ones  of  2007  and 
2009  (Walsh  et  al„  2009;  2011a). 


5.  Conclusion 

One-dimensional  simulations  using  the  ecological  model, 
HABSIM,  indicated  that  late  summer/fall  blooms  of  the  toxic 
dinoflagellate,  Karenia  brevis,  was  initiated  in  response  to  nitrogen 
fixation  by  the  colonial  cyanophyte,  Trichodesmium  erythraeum, 
on  the  WFS,  not  in  response  to  upwelled  or  riverine  sources  of 
nutrients.  Instead,  these  deep  water  and  freshwater  inputs  sup¬ 
ported  a  more  typical  diatom-copepod-clupeid  food  web.  Once, 
K.  brevis  blooms  reached  concentrations  >  4.5  pmol  C  kg-1,  they 
were  capable  of  both  self-shading  and  causing  extensive  mortality 
of  higher  trophic  levels.  The  combination  of  these  decaying  top 
level  predators  and  mortality  of  the  lower  trophic  Trichodesmium 
populations  supplied  additional  non-siliceous  N  and  P  sources  for 
use  by  these  HABs.  This  provided  a  competitive  advantage  for 
further  development  and  maintenance  of  the  slow-growing  K. 
brevis  population  over  the  faster  growing  diatoms  and  micro- 
flagellates.  While  a  monospecific  bloom  of  K.  brevis  at  concentra¬ 
tions  >  5  x  105  cells  1~3  would  take  advantage  of  any  and  all 
nutrient  sources  they  encounter,  the  causal  model’s  specified 
chain  of  events  was  necessary  to  reach  the  observed  cellular 
concentrations  of  red  tides  off  western  Central  Florida  (Fig.  1). 
We  are  thus  ready  to  implement  this  model,  in  conjunction 
with  COMPS  for  development  of  an  operational  model  of  the 
trophodynamics  of  red  tides,  fish  yields,  and  perturbation 
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responses  in  the  eastern  COM  to  climate  change,  overfishing, 
and  oligotrophication. 


6.  Math  formulae 

Diatoms: 

=  T  rt(dP] )-  ^  WiP,  +  dg,  P,  -dy,P,  -dei  Pi  -di/r,  Pi  -d/iPi 
Microfiagellates: 

=  Trt(dP2)  +  dg2P2 -dy2P'2  -  dtg  P2  -  di//2  P2  -  d/2  P2 
Trichodesmium: 

=  Trt(dP3)-  ^(w3P3)+dg3P3-dy3P3-d£3P3-d/3P3 
K.  brevis: 

=  rrt(dP4)-  ^(w4P4)+dg4P4-dy4P4-d£4P4-d/4P4 
Ammonifiers: 

=  Trt(dBa)+dgaBa-dr,aBa~dmaBa 
ot 

Nitrifiers: 

=  Trt(dBn) + dgnBn-denBn—dmB„a 
Dissolved  iron: 

8d^L  =  Trt(dFed)+d[(Fe  :  Qi_4(£i_4Pi_4)-(gi_4Pi_4)+(y,_46yi_y4^ 

+  dppe — dtXpe  +  £  ‘  Oi  _4(S/7i_4/-yi_y4)]  +  dKfe 

Ammonium: 


8d^4  =Trt(dNH4)  +  d(N  :  C) , _4(o, _4P, -4-gi ,2.4Pi  .2.4  +  V 1 -4P1  -4) 
+  d(N  :  C)b(£BagaBa—gnBn)  +  d/lD0N  +  d(N  :  C)i_4£/y//yi_y4 
Nitrate: 

=  Trt(dN03)-d(N  :  C)i_4g,_4Pi-4  +  d£BngnB„ 

Phosphate: 

— ^ —  =  Trt(dP04)+d(P  :  C)i_4(£i_4Pi_4—  gi-4Pi-4  +  yi_4Pi_4) 

+  d(P  :  C)B(8BagaBa  +  £BngnBn—gnBn)  +  df!DOp 
+  d(P  :  Qi-4C/';’//y,  — y4 

Detritus: 

=Trt(dD)-^:(wDD)  +  dqi_4/i_4-dkDD 
Fecal  pellets: 

^  =  Trt(dZ)-^;(wzZ)  +  dAi_47i_4+dyi_5bi_5/-d/fzZ 
Silicate: 

=  Trr(dSi04)  +  d(Si :  C), ffi, P, -g, P,  +A,P,) 

+  dkdiss(Si :  C)i(Z+D) 


DON: 

-  —  =  Trt(dDON)  +  d(N  :  C)i_4(Oi_4/i_4Pi-4+ji_4yi_4Pi-4) 

+  d(N  :  C)b(maBa-\-mnBn—gaBa)-\-d(N  :  C)i_4(kp)D-i-kzZ) 
—dPooN  +  d(N  :  C)^_4kfdF 


DOP: 

yy  =  Tr f(dDOP) d(P  :  C .)  j  4 ( 0 1  4/1  4 P i  4  t ./ 1  4,1  4 P 1  4 ) 

+  d(P  :  C)b(maBa-\-mnBn—gaBa)-\-d(P  :  C)i„4(koD  +  /czZ) 

~dPDOP  +  d(P  :  C)-l_4kfCiF 

DOC: 

<5t^^  =  TrddDOQ  +  dio^l-i^P-i^+ji^y^P-i^) 

+  d(maBa  4-  mnBn  —gaBa )  4-  dfk/jD 4- kzZ)—dl)D0C  4-  d  kjd  F 
Colloidal  iron: 


ddFec 

dt 


=  Trt(dFec)+d(Fe  :  C),^4(Oi_4/i_4Pi_4+ji„4yi_4Pi^4) 

+  d(Fe  :  C):_4(kDD  +  kzZ)-d/3Fe  +  dscFe  +  d(Fe  :  C)1_4kfdF 


DIC: 

dDIC 
et 


=  Trr(dDIC)  +  d(£i_4P,  _4-g,  _4P,  _4  +  y ,  _4P]  _4) 
+  d(£BagaBa  4-  £B„gnB„ -g„B„)  +  d^y,  4fyi _y4 


Fish: 

CZP  =  T  r  t(dZP)  +  dy  1  _4f yl  _y4Pi  _4f>i  -4  -  dkfdZP—  WfZP 


dt 


Diffusive  transport: 

Diatom  growth: 


Zi  =  [/hi 


min(/z//sate<'^/'-<1«),(N03/!kivo3,i)+N03)), 
or  (NFt4/{kNH4ci)+NFl4}),  ( P04/{fcpo4(i)+P04 }), 
(Si04/{kSio4(i)+Si04}) 


g2  =  [P2] 


g3  =  [P3] 


g4  =  1P4] 


Microflagellate  growth: 

min(/z//sate< 1  -WW(2)))xjv03/{kNo3(2) 4- N03 }), 
or  ( NF/4/lkNH4l2)+NH4}),  (P04/lkPo4(2)+P04)) 

Trichodesmium  growth; 

min  (/z//s„te,1^/'-(3)>),(P04/{fcpo4(3,+P04l), 
or  (DOP/ {kDOp^)+ DOP}),  (Fed/{kFed(3)+Fedj) 

K.  brevis  growth: 

min(/z//sate(1-fe/'“(4))),(N03/{klvo3(4)+N03)), 
or  (NH4/{kNHtl4,+NH4))  or  (OON/{kD()Ni4, 4  DON [ ), 

(PO4/ {kpo4(2i  4P04})  or  (DOP/ {koop(4)  +OOP)),(FeB)/ {kFed(4) +Fed] 


Maximal  growth  rate: 

,,  _  11  p0.0633{T— 27) 

7*1-4  —  7*t(l-4)e 

Ammonifier  growth: 


=  [7*0][min(DON/{/<:DoN,a)+DON}),(DOP/{/(:DOp((!)+DOP})] 
Nitrifier  growth: 


Sn  =  [pn][min(NH4/{kNH4(n)+NH4}),(P04/{kpo4(n)+P04))] 
Fish  decay: 
kfd  =  1.97  x  10_1  T 
Iron  scavenging: 
a/e  =  «ma xFed  (Fed/(Fed  +  kaj) 
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Trichodesmium  lysis: 

K  =  [min(P04/{fcpo4(3) +P04}), 
or  (DOP/{kDOp(3)+DOPj ),(Fed/{kFed( 3)  +  Fed] )] 
/3  =  P3/tnl„e~1354f< 
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